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CHAPTER 1 


INTRODUCTION 


The microwave apparent temperature over the ocean measured by a 
radiometer includes contributions from (a) the surface, (b) the intervening atmosphere 
between the surface and the radiometer, and (c) radiation from the atmosphere above 
the surface and any signal source that might be present, which upon reflection 
by the surface falls within the beam of the radiometer antennas* If the atmosphere 
consists of scattering particles such as the precipitating raindrops with size 
comparable to radiation wavelength, this will further complicate the problem. 

A complete understanding of the microwave apparent temperature from the ocean, 
therefore, requires to know the sea surface characteristics and the atmospheric 
conditions that account for the total apparent temperature. 

The apparent temperature contribution from the sea surface, or the microwave 
brightness temperature of the sea surface as usually called, depends on the surface 
temperature and the surface emissivity. The surface emissivity, in general, is 
related to the surface roughness and the percentage sea foam coverage. The sea 
foam coverage especially at high wind speed, say 20 rr)/s, should be properly 
treated. The investigations of sea foam that affects the surface emissivity can be 
found elsewhere (Doppleman, 1970; Hollinger, 1971b; Stogryn, 1971, Porter, 1971). 

It is our purpose here to pursue only the surface emissivity that relates to sur- 
face roughness effect and assume that the sea foam effect is negligible or can be 
properly accounted for when situation requires. With this perspective, the surface 
emissivity depends on the differential scattering coefficient of the sea in an 
integral form defined by Peake (1959). Thus, an assumed scattering model for the 
rough sea also serves as a model for the surface brightness temperature. A satisfactory 
sea model will, therefore, be capable of predicting both the scattering and the 
emission characteristics of the sea surface. This being the case, the roughness 
parameters of the surface which appear in the model may be studies by examining 
either the backscattering coefficient measured by a scatterometer (Bradley, 1971) 
or the microwave apparent temperature measured by a radiometer (Hollinger, 1971a, 
1971b) or by both (Claassen, et al., 1971). In the case of apparent temperature 



measurement in addition , to sea foam effect the atnjiospheric contribution above the 
,sea surface should also be considered. Since atmospheric contribution under clear 
sky or even cloudy sky conditions at low frequencies, say below 15,GHz, is 
relatively small and constant and sea foam effect can also be properly accounted 
for. (Hoi Unger, 1971b, Stogryn, 1.971), both scotterometer and radiometer can be 
employed to detect sea surface roughness . The sea surface roughness in general 
depends on the, wind blowing over the ocean. Thus scotterometer or radiometer data 
can be inferred to wind speed over the ocean . 

Several investigators have reported radiometric observations oyer wind 
driven seas. Hollingerof NRL (1970, 1971a) has conducted radiometric measurements 
at three microwave frequencies, 1.41, 8.36, and 19.34 GHz, from a tower 
located in 60 meters of sea. His measurements for horizontally jpolarized emissions show 
(1) a linear rising trend with wind speed, over all incijdent angle (2) an increasing 
sensitivity to wind speed with frequency, and (3) a small increasing sensitivity to 
wind speed with incident angle. Vertically polarized emission exhibited little 
sensitivity to wind speed except at 19.34 GHz, and large incident angles where a 
declining trend was observed . . In processing his ^data Hoi linger removed the effects 
of foam from the apparent temperature. Atmospheric contributions were also removed 
from the measurements, so only surface geometry affected the measurements since the 
surface water temperature was nearly constant. Nordberget al., (1968, 1969, 1971) 
have also reported the dependence of brightness temperature on wind speed at 19.34 
GHz.* Their observations were conducted from an aircraft over the North Sea and 
the North Atlantic Sea. They report selected data whichjndicated that the 
brightness temperature at nadir increases with wind speed from 7 m/s to 25 itt/s at 
a rate of about 1 ,2^ K/m/s. 

On the other, hand, past theoretical investigations of the microwave brightness 
temperature from the sea surface either assumed a single surface model for computing 
the differential scattering coefficient excluding shadowing and. multiple 
scattering (Stogryri, .1967; Fung and Ulaby, 1970) or including these effects but 
employed one”dimensional model (Lynch and Wagner, 1970) or employed a facet 
reflection model (Sirounian, 1967; Shifrin, 1969). All these investigations 
except Fung and Ulaby's physical optics model failed to predict the observed wind 
dependence at nadir reported by the above measured apparent temperature data. 

Though the physical optics model proposed by Fung and Ulaby (1970) predicted good 
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agreement with experimental results, the parameter in this theory appears as a 
product of the surface parametors which cannot be easily isolated and related to 
wind speed. For practical applications it is necessary to calibrate this parameter 
versus wind speed experimentally.' 

In view of the above deficiency and a closer description of the sea surface 
consists of small ripples, namely capillary waves, riding on large gravity wdves 
(Pierson et al., 1971), it appears reasonable to consider a two scale model (small 
irregularities superimposed upon large undulations) for examining the emission and 
scattering characteristics of the sea. It has already been noted that such a two“scale 
model allowed closer description of the backscattering cross section of the sea as a 
function of the incident angle. Its wind dependence remains to be Investigated. It 
IS the subject matter of this study to further generalize the two“Scale model for 
predicting both the surface brightness temperature and the backscattering cross section 
as a furiction of both wind speed and look angle. 

In the following chapter we shall discuss certain aspects of the sea that have 
significant effects on the differential scattering coefficient.' In particular, we shall 
consider the surface geometry and the electrical properties of sea water. The 
electrical properties in turn are related to salinity, sea water temperature by the 
complex relative dielectric constant. ■ ■ ^ 

In Chapter 3, a detailed development of two “scale non "coherent scattering 
theory extended from Semyonov's paper (1966) for computing the differential scattering 
coefficient Is given. Under the non“coherent assumption the differential scattering 
coefficient can be shown to consist of a sum of two terms: the first term represents 
the main contribution by the large undulations; and the second term the main 
contribution by the small irregularities: The theory assumed Gaussian surface 
height distributions and Gaussian correlation functions for both scales of 
roughness. The differential scattering coefficient is shown to be dependent on the 
rms slope of the large undulations, the standard deviation of the small irregularities 
end the correlation distance of the small Irregularities. The wind dependence for 
the rms slope of the large undulations is treated In accordance with Cox and Munk's 
measured data.' The standard deviation of the small irregularities Is correlated 
to the wind grow capillary wave spectra proposed by Pierson et al'. Good agreement 
with measured characteristics and better agree'merit than the single surface model 
are demonstrated. 
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In Chapter 4, the effect of clouds and rain on the microwave apparent 
temperature over the ocean is presented. Various cloud and rain models investigated 
by meteorologist have been employed to determine the rise in the microwave 
apparent temperature when viewing downward through these model atmospheres. 

A horizontally uniform si-ratified model atmospheres is assumed for both cases with or 
without the scattering term in the solution of radiative transfer equation. The 
attenuation, absorption and scattering coefficients are computed by using Mie 
theory for some representative precipitation rate. With the scattering term Included 
the firiite difference method proposed by Germogenova (1962) is employed for 
solving the radiative transfer problem. Comparisons of the computed apparent tem- 
perature are also made between the two cases with or without the term when the 
sky condition is assumed to be under heavy rain. It is found that the apparent 
temperature contribution by the scattering effect is significant. Hence, for heavy 
rain case, the apparent temperature computation should include the scattering term. 
Emphasis has also been placed in correlating the temperature increase resulting from 
clouds and rain with atmospheric attenuation. Details are given in Section 4.4. 
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I 

CHAPTER 2 


CERTAIN ASPECTS OF THE SEA 


2.1 INTRODUCTION 


A naturally occuring ocean surface is characterized primarily by the salinity, 
the temperature of the sea water, the surface roughness, the foam coverage, the 
white caps,, and the sprays. The salinity and temperature of sea water specify its 
dielectric properties at microwave frequencies in the range 1 to 300 GHz. At 
frequencies below microwaves, say 100 MHz, sea water is a good conductor with 
very high conductivity due to its dissociated salts (Neumann and Pierson, 1966). 

On the other hand the surface roughness, the foam coverage, the white caps, and 
the sprays are developed, in general, by winds blowing over ocean surface. Since 
dielectric properties of sea water and sea surface geometry are responsible for the 
surface emissivity and backscattering cross-section of the sea, it is appropriate to 
have a brief discussion of the dielectric properties and wind roughened surface 
characteristics of sea water. 


2.2 SALINITY, SURFACE TEMPERATURE, AND DIELECTRIC CONSTANTS OF SEA 
WATER 


The sea water is such a complex solution of salts that it is impossible by 
direct chemical analysis to determine its total content of dissolved solids. The sa- 
linity of the sea water has, therefore, been defined as "the total amount of solid 
material in grams contained in one kilogram of sea water when all the carbonate 
has been converted to oxide the bromine and iodine replaced by chlorine, and 
all organic matter completely oxidized. " Since the relative amounts of all the 
major constituents are virtually constants in the different parts of the oceans, the 
salinity of a water sample can be computed by determining the amount of only one 
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of the major components. Chlorine is present in relatively large quantities and can 
be determined easily and accurately; for this reason the determination of the chio” 
rinlty of the sea water has been standardized. The empirical relation between 
salinity and chlorinity is (Neumann and Pierson, 1966) 

Salinity = 0.03 + 1 .805 X chlorinity . (2.1) 

Both salinity and chlorinity are expressed in parts per thousand, or "per mille", for 
which the symbol "o/oo" is used. Table 2.1 gives the major constituents of sea 
water (McLelland, 1965). 


Table 2.1 

PRINCIPAL IONIC CONSTITUENTS OF SEA WATER (34.4 o/oo SALINITY) 


Species 

Amount (g) 
per Kg of 
Solution 

Molecular 
or Atomic 
Weight 

Percentage 
of Total 
Salt 

cations 

Sodium 

10.47 

22.997 

30.4 

Magnesium 

1.28 

24.32 

3.7 

Calcium 

0.41 

40.08 

1.2 

Potassium 

0.38 

39.096 

0.05 

Strontium . 

0.013 

87.63 


Anions 

Chloride 

18.97 

25.459 

55.2 

Sulfate 

2.65 

32.06 

7.7 

Bromide 

0.065 

79.916 

0.2 

Bicarbonate 

0.14 


0.4 

Borste 

0.027 

10.82 

0.08 


From Table 2.1 it can be seen that sodium chloride (NaCI) is the major contributor. 
The specific gravity of a normal solution of NaCl is unity; thus, there are 58.45 grams 
of sodium chloride in 1,000 grams of solution. With this os a basis, the salinity, S, 
may be expressed in terms of normality, N,as (Porter and Wentz III, 1971) 

S = 58.45No/oo (2.2) 
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Table 2.2 1 ists surface salinities during the Northern Hemisphere Summer (Sverdrup, 
et al.„ 1942). 

Sverdrup, et al ., (1942) also furnishes data on the average monthly ocean 
surface temperatures. Table 2.3 lists this information for the months of February 
and August. 


Table 2.2 

SURFACE SALINITIES DURING NORTHERN HEMISPHERE SUMMER 


Ocean Region 

Salinity (o/oo) 


Minimum 

Maximum 

North Atlantic 

34.0 

37.0 

South Atlantic 

34.0 

37.0 

Caribbean 

35.0 

36.5 

North Pacific 

32.0 

35.5 

South Pacific 

34.0 

36.5 

Indian Ocean 

34.0 

36.5 

(to Southern Australia) 




Table 2.3 


SMOOTHED AVERAGE MONTHLY OCEAN SURFACE TEMPERATURES 

Ocean Region 

Temperature (°C) 


February 

August 

North Atlantic 

"1 to 27 

0 to 29 

South Atlantic 

-1 to 27 

-1 to 25 

Caribbean 

26 to 26 

26 to 29 

North Pacific 

-1 to 27 

Oto 28 

South Pacific 

-1 to 28 

-1 to 28 

Indian Ocean 

20 to 28 

15 to 18 


(to Southern Australia) 


On the air-sea boundary three different conditions exist: the stable condition when 
air is warmer than the sea surface water temperature, the unstable condition when 
the air is colder and the neutral condition when air and sea temperatures are the same. 
There are different wind roughned sea surface characteristics under stable, neutral 
and unstable conditions (Cardone, 1969, Figure 2.4). 
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Many invesl’i gators (Lane and Saxton, 1952, Hasted, et al., 1948, 1958, 
Grant, et al., 1957) have made measurements of the dielectric properties of ionic 
solutions of water at various normalities,temperatures and electromagnetic wave- 
lengths. They employed Debye's (1942) formula for comparing theories with 
measurement. The data of Grant (1957) and Hasted (1958) seemed to indicate that 
the Debye expressions for the dielectric constant was in slight error. The data 
could only be explained by a slight spread of relaxation time, as proposed by 
Cole (1941). 

Porter (1971) investigated various parameters that influence the correct ex- 
pression of Debye formula by summarizing all previous results. Following Porter's 
result the complex relative dielectric constant, € r. Is expressed in the form 



4.8 + 


(gs - 4.8 ) [ M- ( -rr X I o -^ ) ] 

14 2 C Va)'" 


(2.3a) 


6 


// 


f - 4.8 ) ( ) 




0.38 


-a- 


1.96 


i8 cr 

■f : (2.3b) 


Where 

€s = 87.8 - 15-3 N - 0.363 T ± 

As = 3.38-0. IIT+ 0.00 l47T^ + 0. 0173 TN-0.52N±^^ 

cr = 5 N + o. I2TN + 0.04 T ± who/„,e-ter 

X = wavelength in cm same as 
= frequency in GH^ 

[v4 = normality 
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T= temperature in 

= 0,1 IS the average deviation for corresponding experimental values 
(>^ = 0.025 
= 0.004 

Table 2.4 shows a set of computed dielectric constants at two frequencies, 
four temperatures and three salinities. 


Table 2.4 

COMPLEX RELATIVE DIELECTRIC CONSTANTS OF SEA WATER 




284 

e' 

288 

292 

e' e" 

296 

€' e" 

1 

GHz 

33°^°° 

52.75 

39.47 

55.06 1 

38.16 

56.66 

36.85 

57.53 

35.72 

9.3 

35 

52.54 

39.60 

54.80 

38.33 

56.34 

37.06 

57.17 

35.99 


37 

52.33 

39.72 

54.53 

38.49 

56.02 

37.28 

56.80 

36.27 


33 

39.93 

39.30 

43.18 

38.93 


38.21 

m 

37.35 

13.9 

35 

39.84 

39.33 

43.04 

38.97 

mm 

38.26 

mm 

37.43 


35 

39.75 

39.35 

42.89 

39.00 

45.42 

38.31 

mm 

37.51 


From Table 2.4 it is seen that the change of sea water temperature affects 
the dielectric constant, while the effect of the change of salinity upon the dielectric 
constant is insignificant at the above frequency range. 

Figure 2.1 illustrates the graphs of the real and the imaginary part of the 
dielectric constant with respect to varying frequencies ranging from 1 to 20 GHz at 
35 o/oo salinity and temperature at 292^K. 

Figure 2.2 illustrates the graphs of the real and the imaginary part of the 
dielectric constant with respect to varying temperature ranging from 276 to 296^K 
at 35o/oo salinity and frequencies at 9.3, 13.9 and 19,34 GHz respectively. 
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2.3 WIND ROUGHENED SEA SURFACE CHARACTERISTICS 


If the ocean were filled with a completely motionless homogeneous water, 
the idea sea surface with a boundary adjacent to the atmosphere would be a level 
surface geologically and a plane conducting surface electrically. The sea surface 
is, however, not of such a simple ideal nature. It varies all the time due to the 
changes of the forces of external and internal pressure and mainly of the wind 
friction stress at the sea surface. 

Any given wind structure is, of course, not uniform over an unlimited 
region of the sea. The length of the path over which a substantially constant 
wind blows is known as the "fetch". At any given fetch the wave amplitude In~ 
creases with the length of time that the wind has been blowing. After a certain 
time a steady state is reached at which the wave height remains constant. It is 
called the "fully developed" sea. The time that it takes to reach this state of 
equilibrium, known as the "minimum duration", increases with the length of the 
fetch and the wind speed. If the wind blows for a time shorter than the minimum 
duration, the wave height does not depend upon the fetch but only upon the time 
and wind velocity. 

In the case of a fully developed sea there Is a unique relationship between 
wind and sea state. Table 2.5, which is an abstract from a compilation due to 
Wilbur Mork of the David Taylor Model Basin, shows some useful relationships. 

Note the wind speed is given in knots in the table (1 knot = 0.5148 n/sec.). 

Waves of the actual ocean are very complex phenomena. A detailed 
description of waves characteristics of ocean can be found elsewhere (Pierson, 1966, 
Pierson, et al, 1971). Since sea surface characteristics responsible for surface 
emission and backscattering cross-section from the large structures can be described 
in good agreement with experiment by using the slope distribution of the sea surface 
measured by Cox and Munk (1952), the wave spectrum analysis of the large structure 
gravity wave is not considered. On the other hand, the description of capillary waves, 
which is also responsible for both surface emission and backscattering characteristics, 
will follow the wind grow wave numbers spectrum proposed by Pierson (1971) 
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Table 2.5 

WIND AND SbA SCALE FOR A FULLY DEVELOPED SEA 


Sea— General | 

Wind 

Sca^* 

Sea 

State 

Dcscn'plion 

Wind 

Force 

(Beau-. 

fort) 

Description 

Range 

(knots) 

Wind 

.Vcloc- 

ily 

(knots) 

i 

Wave Height (Feet) 

Mini- 
mum 
1 Fetch 
(nmi) 

Mini- 
, mum 
Duiaiic- 
(hours/ 

Aver- 

age 

Signifi- 

cam 

Avg 

Highest 

i 

0 

Sea like a mirror 

0 

Calm 

B 

m 

0 

0 

' 0 

... 

' ... 

Ripples with the appearance of scales are formed; 
but without foam crests 

■ 

Light Air 

■ 

■ 

0.05 

0.08 

O.IO 

5 

* 18 mir. 

1 

Small wavelets, still short but more pronounced; 
crests have a glassy appearance, but do not break 

2 

Light 

Breeze 

B 


0.18 

0.29 

0.37 

8 

39 min 

Large wavelets, crests begin to break. Foam of ' 
glassy appearance. Perhaps scattered white 
horses. 1 

■ 

Gentle 

Breeze 

7-10 

8.5 

0.6 

1.0 

1.2 

9.8 

1.7hrs 

2 

10 

0.88 

1.4 

1.8 

10 

2.4 

Small waves, becoming larger; fairly frequent 
White horses 

4 

Moderate 

Breeze 

11-16 

12 

1.4 


2.8 

18 

3.8 

13.5 

1.8 

' 2.9 

B 

24 

' 4.8 

3 

14 

2.0 

3.3 

B 

28 

! 5.2 

16 

2.9 

4.6 

B 

B 

6.6 

.4 

Moderate waves, taking a more pronounced long 
form; many white horses are formed. (Chance of 
some spray.) 

5 

Fresh 

Breeze 

17-21 

18 

3.8 

6.1 

7,8 

B 

8.3. 



6.9 

Bl 

65 

9.2 

5 


5.0 

8.0 

10 

75 

10 
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2.3.1 Slope Distribution Statistics of Sea Surface 

Cox and Munk (1954a) had made surface slope measurement by optical 
reflection method on the Hawaiian sea at wind speeds up to 14 m/sec. Due to the 
method of measurement, the small scale capiljary waves, white caps and sprays 
were ignored in their analysis. To account for such a small scale roughness riding 
on theJarge undulating surface a different approach Is made and it Is presented in 
the next, section. 

' These authors found the slope distribution of the sea surface to be nearly 
Gaussian with small corrections expressed by a few terms in a Gram-Charlier series 
which represents skewness and peakedness. . The principal axes, x and y, were found 
to lie cross wind and parallel to the wind respectively. Cox and Munk expressed the 
slope distribution function in the form. 

PCE*,2,) = (2tt a, (Tj exp (- Yf 7 ^)) ’‘{'“TCaiC?-')’]- 

^ C i ) +■ 2^ Co4 +'3 ) } 

Where X - and 7 = , '^•3 “ — . The mean 

square slopes regardless of direction were found tp depend linearly on the 

wind speed (measured at a height of 41 feet above sea level with stable sea condition). 
The regression lines and correlation coefficijsnt f have been computed by the rnethod 
of least squares ; 

Clean s(!jrface - . 


= ^6.od3 '+ 1:92 X 

1 0“^ vW 

± 0,0.0 2 

. r = 0.9S6 


= 0.000 + 3. 1 6 X 

10 '^ w 

1 0.004 

r = o. 945 

‘ (2.5) 

0-/+ = -O-..00-3 ,+ 5. K2 . X 

lO"^' w 

± 0.00.4 

r ,= o. 9 86 
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Slick surface 

(Tx* * 0.003 + o. 84 X lO’* W i 0.002 

(Jj* = 0.005 + 0.78 ^ W 1 0.002 

0-x% (7,^ » 0.00 8 + 1.5 6 ^ lO'^W 10.004 

J 

Where W is the wind speed in meter per second. The indicated error refers 
to the standard deviation of observed values from the regression lines. Actual values 
of the ratio were found to vary from 1 .0 to 1.8 with a mean value of 1 .34. 
These variations are due to fluctuations in the' wind direction. Values of this ratio 
close to 1 are expected for gusty winds while steady winds lead to higher values. 

The regression lines for the skewness coefficients C 21 and Cqj are given by 

C-, = 0.01 - 0.0086 W ± 0,63 

‘ . (2.7) 

= 0 . 04 - 0.0033 W ± 0,12 

A large scatter in the data is evident from the indicated standard deviations. 

The peakedness coefficients were found to have the values 
C4Q= 0.40- 0.23 

022 = 0.12 - 0.06 ( 2 . 8 ) 

Cq4= 0.23- 0.41 

independent of the wind speed. Cox and Munk state thiert the measured peakedness 
is barely above the limit of experimental error. The observed values may be due to 
either a variability in wind speed or to systematic errors in their experiment. 

This data shows that deviations from a Gaussian distribution are not large. 
Therefore, the sea surface will be assumed to be Gaussian. Because of this assump- 
tion, only;the mean square slope 0 ;*= O’** + given in Equation (2.5) or (2.6) 

is needed in the application that follow this chapter. 

/ 

2.3,2 Empirical Wave Spectrum of Copillary Waves on the Ocean' 

The Capillary waves are presumably caused by very local eddies in the wind. 
They are comparable in size to the radar wavelength. Recent investigation in radar 
return and microwave emission from the ocean surface confirm that these ripples are 
responsible for both backscattering and emission from the sea. . 


r» 0.7 8 
0.7 

r = 0.77 


( 2 . 6 ) 
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Oceanographers have been inferesfed in the longer wavelength parts of the 
spectrum, for these are the one:, that represent sizable energy transfer, and affect 
ships. Published spectra for the sea almost always are based strictly on the gravity 
wave region. Only recently interest has been generated about the capillary waves 
since their growth has been found to be related to the wind speed. 

Measurements of capillary waves on the open sea are extremely difficult 
since they are riding on swells many orders of magnitude larger. Consequently our 
knowledge of capillary waves is bosed on both theory and measurements in a wind 
funnel. Recent work by Pierson et. al . (1971) has combined the data from several 
Investigators so that the variation in the capillary spectrum with wind speed can be 
determined. The empirical wave spectrum of capillary waves, following Pierson's 
(1971) proposed model, may be written in the form 

4 OS X lo”^ 

3 (K) - 


D 


X 


(2.9) 


Where K is the wave number of capillary waves, D is the wind dependent constant. 
The wind dependency of D is carried by the friction velocity, U# defined by (wind 
stress/water density)^'^^, in the form ^ 

D( )= [1.2-47+ 0.0268 U,, + 6.05 X I Uj ] (2.10) 


when is.larger than 12 ^^/sec. 

The results of comparing wind”wave tank data by Toba (1970), Sutherland 
(1967) and Pierson et. al. (1971) are shown in Figure 2.3. The equation shown 
there is a best fit to the observations and indicates an increasing rate of spectral 
growth with increasing friction velocity. Scales are shown both the and 
the wind speed measured at anemometer height of 19.5 m for a neutrally stratified 
atmosphere. For friction velocity below 12 cnr/sec. the. capil lory waves die away 
extremely rapidly and is of no importance for both radar return and surface 
emissivity . 

The relation between the wind velocity measured at z m above the mean sea 
surface and the friction velocity at the sea surface is expressed in the form 
(Pierson 1971) - , ' ' 
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FIG. 2.3 The Wind Dependent Constant, D, of Cap i 1 lary Waves 
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W ( ) = u* /o.4 Log z/%. 


(2.11a) 


z. = -t A 2 8 ^ lo‘^ u/- 4 43 ^ lO* 

(2.11b) 

where both W and U* are in m per sec. 

Cardone (1969) demonstrated that the air-sea temperature conditions 
influence the wind speed measurement at a specific height, and under the same 
atmospheric temperature condition the wind speed varies with respect to the height 
where wind speed is being measured. To make wind speed consistent with real 
surface roughness, namely friction velocity on the surface, every wind speed 
measurement is needed to specify the temperature conditions, namely stable, neutral 
or unstable, and the height of the anemometer. To compare two measured wind 
speeds of different origin a correction factor is needed to bring two situations 
consistent with the same surface roughness conditions and the correction factor 
depends on the anemometer height and temperature conditions. 

Figure 2.4 illustrated the graph of wind speed variation with respect to the 
height of measurement under three atmospheric stratifications. 


2.4 SPECULAR EMISSIVITY OF SEA WATER 


For a plane sea surface, the specular emissivity, related to 

the Fresnel reflection coefficient, R^(6) , in the form 


6*'’ 16) = 1 


1 Rjiej 1- 


= A 


,(2.12) 


where © is the incident angle . The Fresnel reflection coefficient, 
for an air-sea boundary. In turn, is expressed in the form 


R^te) =. 


cos 


6 — J C, — S I rv ^ Q 


cos 


® d c, — s I rt ^ e 


p _ fer COS e - J e 


e, cos 6 -+- ,J e, — e 
where is the complex relative dielectric constant of sea water, 
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(2.13a) 

(2.13b) 
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FIG. 2.5 The Plane Sea Surface Emissivity 


20 


EMiss(V(T,y e^_(o•) 



125 


115 


105 


95 


85 


0 


FIG. 2.6 The Plane Sea Surface Emissiyity and Brightness 
Temperature at Nadir 
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Since fhe complex relative dielectric constant of sea water has been 
discussed in the previous subsection, the specular emissivity can be computed 
by using (2.12) and (2.13). The computed results are shown in Figures 2,5 and 
2,6. Also shown in Figure 2.6 is the specular surface brightness temperature which 
is the product of specular emissivity and surface water temperature. In Figure 2.6 
the specular emissivity at 4.455 GHz is almost temperature independent. Such 
a characteristic may be employed for designing a radiometer as a good sea water 
temperature sensor once the surface roughness effect and the atmospheric attenuation 
are removed or compensated. 



CHAPTER 3 


A THEORY- FOR MICROWAVE EMISSIONS FROM THE SEA 

» • I 

3.1 INTRODUCTION 

The microwave emission charcclerisfic of the sea has been measured by several 
investigators (Nordberg et al., 1969, 1970; Ross et al., 1970; Hollinger 1970, 1971a). 
These investigators have compared their observations with predictions from the 
geomelric optics theory (Stogryn, 1967b) which uses a single^surface model and 
found some but not satisfactory agreement between predictions and measurements. 

The wind dependence of the geometric optics approach was based on measured rms 
sea slope data presented by Cox and Munk (1954). However, the theory failed to 
predict the observed emission characteristics near nadir and fitted only loosely for 
nadir angles of 30 to 70 degrees. The failure of the geometric optics model to 
account for wind dependence at nadir was first reported by Nordberg, et al . (1969), 
and was verified by Hollinger (1971a). 

In view of the above deficiencies, an investigation is necessary to seek a 
more adequate model for microwave emissions from the sea. The emphasis in this 
investigation is oriented towards using a composite-surface model which better reflects 
fhe roughness characteristic of the sea. Since several lengthy numerical integrations 
are required to yield the emissivity, the more adequate model must not be so compli- 
cated that it makes numerical calculations prohibitive. With this perspective, a non- 
coherent scattering theory of the type described by Semyonov (1966) is extended to 
yield the bistatic scattering coefficient. Since an acceptable scattering coefficient 
for predicting the microwave emission characteristics must also be acceptable for 
predicting backscattering, backscattering is examined to provide a cross check on 
the model . 

To compare with experimental observations, an Isotropic surface characteristic, 
although it is not realistic for the ocean surface, is assumed. A justification for 
this assumption is based on the observed directional insensitivity of emission from 
the sea (Hollinger, 1971b). The two-scale rough surface model is also assumed to 
hove Gaussian surface height distribution and Gaussian surface correlation for both 
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scales. The dielecfric constant needed rri the calculations is based on the data 
reported in Chapter II. ♦;* 

The wind dependence of the surface parameters in the composite model is 
introduced In accordance with r ms slope data measured by Cox and Munk Tor the 
large undulations and Sutherlands (1967) results for the small irregulbrities. Details 
for the theory and the choice of surface parameters are given in later sections. 
Comparisons of the computed brightness temperature and backscattering characteristics 
at two different wind speeds are made with both measured data and the predictions 
of a single surface model . 

. . . . . ’ ‘y 

3.2 SURFACE BRIGHTNESS TEMPERATURE THEORY 
3.2.1 Formulation of the Problem 

^ The basic theory of apparent surface temperature was developed by Peake 
(1959). The study that fojiows concerns only the contribution to the. apparent 
surface temperature due to thermal emission from the sea surface; contributions due 
to sky radiation reflected from the sea are not considered. The relationship among 
the surface emissivity, the surface temperature and the apparent surface temperature 
(or the surface brightness temperature in this case) is as follows; 

where is the brightness temperature; the emissivity; Tj the surface 

temperature; ft. Is the horizontal polarization and v is the vertical polarization. 
Note that the azimuthal angle 4> needed together with the nadir angle © to 
specify a direction has been chosen to be zero without loss of generality. 

The connection between the emissivity and the differential scattering coef- 
ficient of the surface, ) is 

/27T r Va 

''o ''p ^ 

where 0^ , and % are angles' defining the direction of scattering corresponding 
to a wave incident at an angle 0 . 
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The basic formulation of the problem indicated above shows that the different 
fial scattering coefficient is the quantity that defines the angular characteristics of 
the brightness temperature of a given surface. Consequently, different brightness 
temperature theories are also distinguished by the different models assumed for the 
differential scattering coefficient. 

Under the non-coherent assumption ) can be shown (Semyonov, 

1966) to consist of two terms, 

‘ ( 0 . 9s, ) = V ^ ^ 

where <ti)clenotes the main contribution by the large undulations; and 

< Y^'C9, denotes the differential scattering coefficient of the small irregularities 
averaged over the distribution of the surface normals of the large undulations. 
Detailed derivations for V.“C©/ and < y.'c®/ ®». 4*s are given in the follow- 

ing two sections . 

Since the backscattering cross-section is a special case of the differential 
scattering coefficients, it can be obtained from the known differential scattering 
coefficients, (see Figure 3.1) i,e. 


a" (8) = cos e • Y^(6, 


or equivalently 


(6)= <^Bjo (0) + ISj) 


(3,3o) 




(3.3b) 


6 (y/(0.es.4>s)> 


(3.3c) 


Detailed derivations of <S^AQ)ore given in the section on backscattering cross 


sections , 
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3.2.2 Derivation of y/i®, ‘t’s ) 

To derive ®». Is) we begin with the vector-scattered field due to o 

plane wave incident on an undulating surface to which the tangent plane approxi- 
mation is assumed applicable. Such a field expression is, in general, rather 
complicated. To simplify results, the stationary phase technique will be employed. 

An expression for )< S, “ft ) not indicating explicitly the effect of the small 
irregularities will be derived first. This expression is then clarified to reflect the 
small structure effects by computing the explicit forms of the modified Fresnel 
reflection coefficients. As pointed out by Semyonov, such a computation may be 
performed for the more general finitely conducting surface in accordance with the 
work of Rice (1951). 

(i) The Scattered Field 

The far zone scattered field in the direction 3k* clue to a plane wave polarized 
along ^ impinging on an undulating surface can be shown to be (Fung, 

1967) 

Ea.s = )CzkJ<.:b ) + (» ) t] 

-»■ 2^2 X [c ( a • ^ Xzi 5<i ) -V C \ - <R<,» (5, • c a ' a . ) t 

MM 

• E. exp[-a(Us,- ao-r] ds. „ 



J U)t ~"J ^ ^ / 

where a time factor of the forme has been suppressed; X = -;ke /i4n/sj, R is 
the distance from the origin to the field point, k is the wave number; <Rtr'> ce 
the modified Fresnel reflection for horizontal and vertical polarization respectively; 

Eo is the magnitude of incident field, and is the normal to the surface . 

The set of orthogonal unit vectors (jii^i^d,) serving as the local coordinates 
for evaluating the local field on surface is illustrated in Figure 3.2. The unit vectors, 
t and d relate to 2, and as follows 


t = C *i. X n)/ I n ( 

d ~ 21« X t 


(3.5) 
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Figure 3.1 The sefs of orthogonal unit vector , if;) and 

i :: ( n,/ I^t) for expressing the polarization states of 

th^ incident and the scattered field. 


Z 



d 


Figure 3.2 A sef of orfhogonal unit vectors , i, , ^ ) serving 
os the local coordinates for evaluating the local field 
on the surface . 



where may be wriffen In terms of the partial derlvaties, of Z(x, ^ ) 

along x and y axes as follows 

) 

n = C-i ?x - i i 3 ^ ^ . ' 

( ^ JL , tariff coorcl^^a^e vecfors. 

To simplify (3.4) by fhe stationary phase approximation, let ^ 
and ^ be the vector components of f . Then the phase factor in 

(3.4) is ' ' "" 


I'l = 


(3.6) 


The stationary phase assumption requires that 
It follows that 

2x = -?./?, ,.v - 

^3 t)/?f - (3.7) 

2 i = l/lll 


The significant result indicated by (3.7) is that all surface slopes in (3.4) 
can be written in terms of the incident and scattered propagation vectors. Conse- 
quently, the integrand in (3.4) except for the phase factor, evpr- j ^ X.' ^ ^ 
be moved outside of the integral sign; i.e., 

= KKZ}:z X ^ j 

. Ao ,, 

(3.8) 


with 


A = (* 

A = 0 - <R:>> ) ( ^ • i ) ( n • n. ) t 
^ +<R,>)( A- c()( n. X t ) 

P = ( 1 -<R^>)(^- X) Cn ■ n, ) t 


(3.9) 
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The local coordinate vectors, two other vector products In (3,9) and the dif- 
ferential surface element in (3.8) can all be written in terms of the propagation 
vectors, 

; i X nO /l a, x.n:, 1 ■ . ■ 

^ i a. X Til I 

■ 2 1 X I 


■n. • Hi = - ' 

^ ^ 2 

.. j ■ ■ ' 

dS = Cl tlj - n.,\ ) dxd^ 


(3,10) 

I 


To express the polarization states of both the incident and the scattered 
field, it is convenient to introduce a set of orthogonal unit vectors ( ' ) 

for the incident field and another set ( ;jTt ^ for the scattered field (see 

Figure 3.1), In view of Figure 3.1, explicit expressions for these unit vectors may 
be written as 

"“Z!:' “ ^ J;. Cos ^ ^ 

^ <^Ps G - C, -h 3 in. 9 . ^ 

'(3.11) 

i,= X ■ . ^ ^ 

= -sin. 04 COS 4’s sm siVv 4>s -4- Cos. S, 

iti = - ^Sds. cos 't’s - cos 0^ Sin <k 'i S,V &. ^ 

•cL ^ 

si^ 4s 3 (3 »12) 


For horizontally polarized incident wave, (£^= 1,), the polarized and cross 
polarized scattered fields are given respectively by 
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= t„S 




)j ]| ^xp[--a^(ii~st)-l]^^ 

(3.13) 


= EcK[ 


u, • n j X 


[( £L{.~ 




2v)l]j 


a,-tt.)-z]ds 


where 


A,^ - \^\i\ -V <Rfc>)( 

'w ^ 

l^Ui - <Rv>Hi,-:riO(:i'>' ^i») 

C _-\ 4 ^ 1 (\ -+-<R.»(;ir 2 lOC S;.+J^O 

— «^H 

D^= l-fK 1 - - 2 ;Od' 


Ifl = \ :^l-v^^\ /2\ u.'X 

V 

Similarly, for a verHcally polarize'd incident wave ( & = polarized 

and depolarized scattered fields are 


^vk -“ 

= E. St i .- + 5 . X C£v + SvM ) • 


[ „pt j^tni-ao-x^ ds 

■*A. 


= £* ' I 


(3.14) 


V,s 




:k* x( c„ ^ 2 ^)]] 


.[ e;.p [}{c(Zlt-3:i)-l] 

-'A. 
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v/here 


l-fl ( ' -V- <R«»(Si- 
B„ = lfl (1 -<R^>)Cv, 


=-lfl ( \ <i?v>)( :ir. • 5^)C2:. 2^t) 

'Vvv'^ 

D„=-WO - <R«>Ki'-;iO(2'XaO 

<Vv 

These field expressions can be simplified furfher by using the vector 
identifies 

X ( 21. X N ) = U,-M - -N 

ir.- X c M -V )i. ) = - i-, - M - or,.N 

Thus, 

C, ( -<R^> be + <R„> d-f ) I 
Ei„= C (<Rt> bf + <R^> dc)j 
^vi.^ CA- di 
Evv= C, C <K^> dc -h <Rv>bf ) I 


where 







®s <^°S — s'in. 0 

^3 ~ ® ©S S'l rt 

= Coi 6 -V cos 0^ 
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b — ( V, ' ^ 9s Cos 0 s»vt 6s cos 

e . = ( Ui • 2:t ^ ~ ^ ^ s\a 03 ' 

d =5 ( ;^, • == s\w 63 s'\A 4s 

“f =■ ( ^1 * ^’"^.9 si A 4^s 


^1 1 iJ -2 ^'V:\ I V^2 = ' 1 -h Cos 6 Cos 0s ~ Os cos'^'s' 


Nofe that by expressing in ^erms of fhe orffiogonal sef ( v;, ^ ) 

if can be shown fhat 1 21t ^ Ji' 1 * = b* ■+ d ^ . Similarly, by expressing , in 
ferms of fhe orfhogonal sef ( iJ3u, ) if follows fhaf 1 S' x Si 1*= c* -v -f *. 

(ii) The Differenfial Scaffering Coefficienfs 

The differenfial scaffering coefficienfs relafed fo fhe scaffered fields computed 
in fhe previous section (Peake, 1959 and Storgryn, 1967a) are of the form 



(0, ds.^s) = 


4-Tt K Bjc £).- > 
^ A, c,os e 


(3.16) 


where fhe symbol denotes complex conjugate, and > the ensemble average. 

The subscripts j, i denote the polarization states of the incident and the scattered 
fields respectively. is the Illuminated area. 

Upon substituting (3.15) into (3.16) it follows that 


V/ ( 6^ 6s, ) = ya ( 5/ ®s; J + y<ui ) 

/ |T|^S 

~ TiA,,cos& f/ i + d" J 


(3.17a) 


( 9, 6s, 'fs ) = y^i ( 0, fs ) -i" Vv ( ®/ “Ps ) • • 

r , 2 ’ (3.17b) 

_ C a,^ f ( < III' > 

irA^cosb'l^ [ b' + d^ 

• \ 
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For an isotropically rough surface with Gaussian height distribution, <l'X(^ > 
simpi ifies to 




(3.18) 


where ^is the zero order Bessel function of the first kind; and ) are 
the variance and the autocorrelation coefficient of the surface respectively* 

An approximate solution for (3.18) consistent with the stationary phase 
approximation is 




(3.19) 


where m. = C <5" l f 'co)l ] ^\s the rms slope of the surface. 

3.2.3 Derivation of < 

\ I ■ 

The scaffered field due only fo the small irregularities has been derived by 
many investigators (Rice, 1951, Valenzuela, 1967, and Barrick and Peake, 1967). 

As was indicated in the previous section, when the scattered field expression is 
known, the scattering coefficient can be computed. To account for the interaction 
between the small Irregularities and the large undulations, the expression for the scat- 
tering coefficient of the small irregularities is averaged with respect to the slope 
distribution of the large undulations. The resulting expression is the desired 

scattering coefficient, < Y-\ 6, 6», 

3 , 

(i) Differential Scattering Coefficients 

The far zone scattered field of i. -polarization along the direction defined by 
the angles 6^ and 4s^ due to a plane wave of ^ -polarization with unit 
amplitude impinging on an irregular surface along the direction defined by 

the angles 0' and has been derived by using the method of small perturbation 

(Barrick and Peake, 1967). The ensemble average of the magnitude square of the 
scattered field, O ( S(, S/ j fs ) 1^ can be shown to be (Barrick and 

Peake, 1967) 


<1 I W c p, ) 


(3.20) 
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where 


- incident and scattered polarization respective, either horizontal 
or vertical polarizations, 

*^0 = illuminated area, 

R = distance from the field point to the surface, and 
M., - C£r - I ) cos cp') 

[ Cos e% Ce, - [ cos -+■ ( e' J 

_ («-.-!> sin - ‘i’') C^r - SiK^ 

[ WS eV C6, - s.'w' e') J [ cos Q' (e,~ 9j") ] (3.21) 

|v< _ ~ - I ) Si'"- C«t’s - (€. - 6s ^ ^ 

[e, cos ] [ cos &s -\- ( e, - siV’Sj ^ ] 

1 ^ _ - Ce. - 1 1 fc°s c t') C€r- c^r- )^ - e, s.^e,' siVt e'J 

[e, Cos e'-i- (e, - sin .2 9 ' V ''2 3 [e, cos 65 -+ ce,- e^; ) ] 


in which is the complex relative dielectric constant, 

^ = surface roughness spectral density, 

P = ^ [ siH. 64 00 s ( 4^'- 4«') — si'w, 9'') 

% ~ s>K 9j^ siK c 4>s — 3 

Note that the surface roughness spectral density W ( p, | ), is related to Its 
correlation function R(^,^)by 

-u» 

For an isotropically rough surface, it reduces to 

RLr)==:a-J WU)T,Ctr)tdt 

o 

and 

Wlt)= Af^p^r) 7 „Ur)rdr 

't> 

For a Gaussian spectral density, it is expressed in the form 

. (- (3.22) 
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where 




t = V) 

^ y i = sfandard deviafion and correla^ion length of the small 


irregularities respectively. 

From (3.16) and the relation 

( »/ 4 ; e,: 4,-j = y.^ i e/, <fi ) 

we get 

j = i.,v 

(3.23) 

(ii) Averaging Procedure 

To account for the tilting effect of the small irregularities by the large 
undulations it is necessary to average ( 0 / 6s ^ ) with respect to 

0 

the slope distribution of the large undulations (Semyonov, 1966), That Is 

<'n\6,6s,%)y =j f 

i ^3^24) 


where 6 is the incident angle and 9s and are the scattering angles. 

To evaluate the above Integral, connecting relations between the surface 
slopes and the local angles, 6', <Pj 9/^ are needed. To find 

these relations let us first express and in terms of the azimuth and the 

elevation angles, 4 ^ and 9^ / which represent the tilting effect: 

7 = cos <}> tan 6n 

= s/'« (3.25) 

1 


From (3.25) It follows that 



3ec ^ 0^ SIM. 0^ d dn '^n 

35 


d Z^ — 


(3.26) 



The local angles (o' ^ ) can now be relafed to the azimuth and the 

elevation angles and ©a by connecting relations derived below. 

In Figure 3.3, the two sets of coordinates {x ^ ^ ^ ^ ) and { , }^ ) 

are related in terms of the angles 9^ and 4^^ as follows: 


X 

tl * 


? 



cos 9^ cx?s 

~ 0^^ Cos 


cos 0K 

cos 4vt 

— s iVt Qrt. s I vv 4^ 


o 

cos 




(3,27) 


Hence, for a scatfered field poin^, P, locafed at a distance, R , from the origin, 
the coordinates of P may be expressed either in terms of the angles 9^ and 4>/ 


lies 

05 

and 

r 


x'=s 

R 

s\n < 

5s cos 4>s 



R 


6s si'^x 4’s 

(3.28a) 

5' = 

, R 

cos 
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X =» 

R 


Ss *^5 


a = 

R 

sr^t 

®S S**'- 

(3.28b) 


T?. Cos 


Substituting equation (3.28) into (3,27) we obtain the connecting equations for 
the sets of angles as follows: 


"*Sin. 9 ^' cos. ' 


fcos c°se« 0H 1 


" Gj Cos 4>j 

Sivv 9^ ^‘iv, 4^ 

= 



SiH Oj S‘«n 4 *^ 



Si'k 6^ cos ^ ^ i*' ®« ®*« J 


0s ^ 


(3.29) 


If we take the angles with a prime to be local scattering angles, we see 
from (3.29)' that the local scattering angles may be expressed in terms of the scat- 
tering angles 85 , and the tilting angles , i.e. 


cos fij = cos cos 0s — SlK-. Sn si^^ Ss Co S ( 4>s - ) 

si'~ 4 ^' = siK Ss s'.« ( 4>s - 4*.> ) C I — 




(3.30) 
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In Q similar fashion, fhe local inciclen^ angles can also be expressed In ferms of 

fhe tilting angles 0^, 4^, and the incident angle, 0 . 

{ 

cos 8 =, Cos 0» cos & -f- s.v 6^ stK 6 cos 4 'k 

sfh. - siv, G- si'n - cos* e ' (3.31) 


When the connecting relations between the local angles and the surface 
slopes are known, (3.24) can be evaluated by assuming a form of P( ) .■ 

For a Gaussian surface slope distribution it may be represented as 


P(2x 




I 




ex 


p ) / 2 j 


(3.32) 


or equivalently 

exp (- J (3.33) 


where m, the rms slope, is assigned according to Cox and Munk's slick sea data. 

Since m.* is usually sufficiently small for the sea and since y.’ 6s 4i' ) 

is Insensitive to 0^ for small values of Qn, the integration with respect to 0,^ 
given in (3.24) can be evaluted by the -method of steepest decent, the result being 
expressed as ' ' 


<Y.'(6, )> 






-Yz 


\ 4>: 64, 4-,') 




v/(&; 


Cm; 


d 


(3.34) 


3.2.4 Modified Fresnel Reflection Coefficients 

As was mentioned earlier, the Fresnel reflection coefficient should be 
modified to account for the presence of the small Irregularities (Rice, 1951, 
Semyonov, 1966); The method for computing these coefficients has been discussed 
by R ice for horizontally polarized waves. Following Rice's approach Valenzuela 
(T970) derived the modified reflection coefficient for the vertically polarized waves. 
However, there appears to be an error in his results, since they differ from Rice's 
original work by a cosine factor in one of the terms. The corrected modified 
Fresnel reflection coefficients calculated by Rice's method are 
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where 


K (®; = cos © ~ (gj, - sih.^® ^ ^ 
Cos ® -f- (e^- 


P /@) g, Cos® - fe, - s/;.^( 
^ Cos s> -h Cei-- 


(3.36) 


Cos © = 


f r - _ . . iiT 

C I + cos © cos 05 — Sl« 5 s/m 0s C®S ^ 


in which Sp is. the complex relative dielectric constant and where 


H. 


e « C u.^ + v-')] 

b Ti.’’- C«‘ + u") J ^ 


(3.37) 


For the purpose of this paper the surface roughness spectrum in the above 
formula is chosen to be 

W Cu- -fes.K ®, vj = ^ exp|^--if (u’ + v^- liusin ®+ 


(3.38) 


3.2.5 Backscattering Cross~Sections 

Substituting 65 = 6 and 4s = into (3.36) we obtain 


R<.Co ) => = 


( - U. 1 




1 ^ 


(3.39) 
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From (3.17), (3.19) and (3.3b), it follows that 


^Buo(e)= 


(3.40a) 


(3.40b) 


To find < )>■, note the following two points; 

(1) Substituting %= 6 and %= 'fl into (3.30) we obtain 


COS ^ a=. Cos 9 =; cos cos 0 Sj'k 0^,COS Sivv 0 
4>5 = - S1.V e^>' =, s,'w 


(3.41) 


or 


8 s' = e' 


= 4»' -V- 7r 


(2) Since 85 = 6 and 4^ = 7r v; imply 4»s = 4 >Vtt and O 5 = 6 ^ , 

from (3.34) it follows that 
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(3.42> 


Thus frorri (3.3c) and (3.42), 


<cT3 .,IG)> = 


2 K 






(3.43) 
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where 


V = 0 ' 

= 4-1'^ 0-," j5* cos'^ &'j Rf^Ce'^l^ exp (.- ■!.-.n' &' ) 


(3.44) 


and 


BUI 


(e-) = c... e' v‘ (.e; e,', 4>.')| 


s; = 9 


— 4-fe'^ cr* Ji* Co^ e' U.'^>--I n &' t- e, 3 1 ^ exp(-vrsi«* 

. 1 r^r eos^.G/'-t- ^ ‘ 

. . :• (3.45) 


The complefe backscattering cross-section Is, of course, given by (3.3a) that is 
the sum of (3.40) and (3.43). ' 


3.3 SELECTION OF PARAMETERS 

The surface parameters that appear in the above theory are the rms slope 

of the large structures m , the standard deviation of the, small structures , O’, , 

and the correlation' length of the small structures i . This scattering model can be 

adapted to predict sea returns by notir^ the rms slope can be based, on , 

measurements by Cox and Munk (1952) and that the assumed Gaussian spectrum for 

the small irregularities can approximate the high frequency port of the sea spectrum 

BK where the Bragg scatter condition applies, i.e., K= 2k sin 0 . In view 

of the requirements of the composite surface theory it Is reasonable to choose the 

oil slick sea measurements by Cox and Munk, slnce the small/ irregularities hove 

been suppressed. The value of m Is thus assigned. The value of kl is assigned 

-4 

so that the correct angular behavoir of the Gaussian spectrum approximated BK 
well over the angular range, '30^^ ^ ^ 70^ i.e., BK ^ is approximated by 
4^/ir where K = 2k sin 6 , the Bragg scatter condition. This 

approximation is achieved by noting that when kZ = 2 similar behavoir is realized 
(see Figure 3.4), The factor 35.3 appears in the Gaussian approximation to 
bring the levels into agreement at 0 = 60 degrees. The value of B must yet be 
assigned. 


41 





Oceanographic inves^iga^ions indicafe values of B in l^he range from 4.6 x 
10 ^ fo 3.26 X 10 ^ (Cox and Munk, 1954, Pierson, 1970, Phillips, 1966), This 
implies fhaf should lie in the range from 0.084 ^o 0.24 when BK ^ is equafed 
fo ('K. J?%_} af 0=6Odegrees. These values of k <;r^ are cons is fen f 

wifh fbe assumpfions of fhe small perfurbafion fheory, an encouraging result. The 
recent reports by Sutherland (1967) and Pierson et al , (1971) indicate that B is a 
function of the wind speed. Thus, the surface parameter must also be a 

function of the wind speed. It Is noted that the horizontally polarized emission 
characteristic for nadir angles from zero to thirty degrees is very sensitive to kcr, 
and hence the parameter k o", can be estimated by fitting the predicted emission 
characteristics to the measured data. (It appears that the wind sensitivity of B may be 
assigned by this technique.) 

With the surface parameters established in the manner described above, both 
the emission and backscatter characteristics may be computed and compared with 
reported measurements. 


3.4 COMPARISON WITH EXPERIMENTS 


The parameter ka^ is estimated from horizontally polarized emission character” 
istics at 8.36 GHz associated with two distinct wind speeds. The emission characteristics 
are based on an average of several of Hollinger's experiment runs ‘under similar 
wind stress conditions (Cardone, 1969). The data reported by Hollinger are in the 
form of surface brightness temperature. Effects due to foam and the reflected sky 
radiation have been removed (Hollinger, 1971a). The vertically polarized emission 
characteristic is computed from the estimated ko^ . These results are shown in the 
graphs of Figures 3.5 and 3.6. The dielectric constant is based on data reported by 
Saxton and Lane (1952). 

Comparison of the predictions of this emission model indicates a significantly 
Improved agreement over that predicted by a single-surface model. Better level and 
trend agreement Is evident for both horizontally and vertically polarized emissions. 
Sensitivity to wind speed is evident at nadir which is not noted with the single-surface 
model. The sensitivity at nadir is carried by the modified Fresnel coefficient (see 
Equation (3.35)). 


43 



( 


'\- 


The adequacy of fhe composlfersurfacc. theory is further demonstrafed when 
the predicted backscatter characteristics are compared with measured characteristics. 

• Data at 8.91 GHz reported by Dale))', et al ., (1971) at similar wind stress conditions 
were chosen qs a basis for comparision. The dielectric constant is changed to 
reflect the influence of the slightly different, frequency . The comparison of predicted 
and measured characteristics are shown in Figure 3.7 to 3.10. These results indicate 
reasonable^dgreement with measurements and improved agreement over the predictions 
of the simple geometric optics approach (single-surface model). It is noted that the 
best agreement with measurements occurs primarily at larger angles and for vertical 
polarization. There is some uncertainty in the accuracy of the measurements near 
nadir (Daley et al., 1971) so lack of agreement may be anticipated there. The 
discrepancy at large angles for horizontally polarized cross sections may be attribu*" 
table to receiver noise at these small cross sections. This statement is, however, 
speculative. " ^ 

The level of NRL data which are based on the statistical median had to be 
raised by 6 dB to realize the agreement. Valenzuela et al., (1971) indicated that the 
overage cross section was about 4.6 dB above the median based on exponential 
statistics assumed for the returns. As a consequence, 1 .4 dB remains unaccounted 
for. Perhaps the remaining 1 .4 dB may be partially associated with the biases 
disclosed by Claassen and Fung (1972), 
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0 10 20 30 40 50 60 70 


I NCI DENT ANGLE 0 (DEGREES) 

Figure 3.5 Comparison of surface brightness temperature between the measured 
data and theoretical calculations using the composite~surface theory 
and the single-surface theory (after Stogryn). 
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I NCI DENT ANGLE 9 (DEGREES) 

Figure 3*6 Comparison of surface brightness temperature between 

the measured data and theoretical calculations using the 
compos ite*“surface theory and the slngle^surface theory 
(after Stogryn). ‘ 
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INCI DENT ANGLE e (DEGREES) 

Figure 3.7 Comparison of backscaffering cross“secHon between the 
measured data and theoretical calculations using the 
composite-surface theory and the single-surface theory 
(after Stogryn). 
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BACKSCATTERING CROSS SECTION o°o(db) 


8,91 GHZ. HORZ. POLARIZATION 



0 - 10 20 30 , 40 50 60 70 80 90 

INCI DENT ANGLE 0 (DEGREES) 

Figure 3.8 Comparison of backscattering cross-secfion bel-ween |-he 
measured data and theoretical calculations using the 
composite-surface theory and the single-surface theory 
(after Stogryn). 
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BACKSCATTERING CROSS-SECTION 



0 10 20 30 40 50 60 70 80 90 

INCIDENT ANGLE 0(DEGREES) 


Figure 3.9 Comparison of backscaH’ering cross~section belweeh the 
measured data and theoretical calculations using the 
composite-surface theory and the single-surface theory 
(after Stogryn).^^ 




Figure 3.10 Comparison of backscattering cross-section between the 
measured data and theoretical calculations using 
the composite-surface theory and tlie single-surface 
theory (after Stogryn). 
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CHAPTER 4 


THE EFFECTS OF CLOUDS AND RAIN ON THE 
MICROWAVE APPARENT TEMPERATURE 


4.1 INTRODUCTION 


High winds over the ocean are usually accompanied by clouds and rain 
(Claassen, et al., 1971). These meteorological conditions influence radiometric 
observations over the ocean primarily because of atmospheric absorption (emission). 

In the particular case of heavy rain fall, radiometric measurements are affected 
also be scattering due to rain drops. It is the aim of this chapter to discuss these 
meteorological effects on the measured or computed apparent temperature as observed 
by a radiometer. 

^ In clear sky the contributions to the apparent temperature are mainly due to 

the absorption of the molecular oxygen and water vapor present in the atmosphere. 
The absorption ot these constituents can be calculated by using the well developed 
formulas (Meeks and Lilley, 1963; Staelln, 1966) and the model atmosphere which 
describes the temperature, the pressure and the moisture profiles (Valleys, 1965). 
Hence, in clear sky the apparent temperature due to the^ atmosphere is reasonably 
well defined. 

In cloudy sky, the liquid water content of the cloud absorbs and hence emits 
microwave radiation in addition to the effects of the molecular oxygen and water 
vapor already present under clear sky conditions. Furthermore, different kinds of 
cloud with various spatial distributions may be present in the atmosphere. Con- 
sequently to investigate the cloud effect on the measured apparent temperature is a 
much more complicated problem than the clear sky case. 

In rainy sky the raindrops not only can absorb but also scatter microwave 
radiation in all directions. If the size of the rain-drop is much smaller than the 
incident wavelength, it has been shown (Kerr, 1951) that the scattered energy is 
negligible as compared to the absorbed energy. The terms representing scattering 
etfects in the radiative transfer equation can therefore be ignored. However, when 
the size ot raindrops Is comparable to the Incident wavelength, the scattered energy 
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becomes significanf and the scattering effect has to be included. In the later case 
a mathematically Involved formulation is needed to solve the problem. 

The apparent temperature contribution to the surface brightness temperature 
by clouds and rain have been reported by Singer and Williams (1968) at 15.89 GHz. 
Additional observations at 19.45 GHz were reported by Conaway (1969) and Kreiss 
(1969). On the other hand, theoretical investigations of cloud and rain effects, 
except for the efforts of Kreiss (1969) and Porter (1971) have. largely been limited to 
cases in which only the sky temperature is computed (Stogryn, 1964; Shifrin, 1969). 

Other radiometric measurements have documented the effects of clouds and 
rain on the sky temperature^ i.e. looking upward. For example, Haroules and Brown 
(1969) in a comprehensive set of measurements have measured sky temperature and 
inferred atmospheric attenuation. The attenuations and sky temperatures of various 
cloud and rain conditions shown in Table 4.1 are an excerpt of their data. Other 
extensive measurements of this type have been reported by Strickland (1970), 

Crane (1971), and Ippolito (1971) in conjunction with the ATS-V communication 
satellite experiments. From these results the apparent temperature when looking 
downward may be qualitatively described. 

Radiometric measurements conducted over clouds and rain will experience 
an attenuated surface brightness temperature plus a cloud or rain contribution in 
excess of the clear sky temperature. When the attenuation due to clouds and rain is 
sufficiently large it is possible for the radiometer to lose contact with the surface. The 
greater sensitivity of the rodlometer to clouds and rain can be employed to correlate 
the increased apparent temperature with attenuation when clouds or rain intervene. 

This attenuation measured from the excess apparent temperature, in turn, can be 
used to calibrate scatterometer measurements and to determine when scatterometer 
data should be discarded due to the presence of heavy rain. Although the relation 
between atmospheric attenuation and apparent temperature has been studied for 
upward looking radiometers, and measurements have shown that the expected 
temperature enhancements do Indeed occur for downward looking radiometers, this 
is believed' to be the first study aimed specifically at determining the relationship 
between brightness temperature and attenuation for the downward looking radiometer. 
Apparent temperatures for various meteorological conditions are computed by 
employing the models postulated by meteorologists. Comparison of the apparent 
temperature with attenuation has justified the general idea of the use of radiometer 
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TABLE 4. I 

ATTENUATION AND BRIGHTNESS TEMPERATURE OF CLOUDS AND RAIN 


Cloud Type 

One-Way Attenuation 

^sky 

8 GHz 
dB 

15 GHz 
dB 

8 GHz 
, °K 

15 GHz 

High Level 

.02-. 05 

.03-. 07 

7 

15 

Medium Level 

.03-. 12 

CO 

1* 

CO 

o 

10 

21 

Low Level 

.10-. 25 

1 

o 

14 

30 

Light Rain 

.15-. 25 

.20-. 60 

16 

38 

Medium Rain 

o 

00 

3.2 

20 

43 

Heavy Rain 

1.5 

5.4 

70 

140 

j 
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to calibrate the scatterometer . The surface roughness model needed In the 
computation is based on the composite surface model developed in Chapter III, 


4.2 THEORETICAL APPROACH 

4.2,1 The Radiative Transfer Equation 

To study the radiative transfer problem ot the planets, astronomers have 
formulated the equation of radiative transfer. By solving this equation the intensity 
of radiation from the planets under specific conditions can be obtained, A down" 
ward looking radiometer at an altitude above the ground surface receives the 
intensity of radiation from the ground surface as well os the intervening atmosphere. 
Hence the radiative transfer equation can also be applied to compute the intensity 
of radiation received by the radiometer. The radiative transfer equation for radio- 
metric temperature computation has already been formulated (Stogryn, 1964; 

Shifrin, 1969). In radiometric temperature computation, the apparent temperature 
at any point Is defined os the temperature of an equivalent blackbody which 
emits the same intensity as the existing radiation intensity at any point. Thus, at 
microwave frequency, the apparent temperature e, <P ) is related to the 

radiation intensity in the form 

T.jCj, 0, (4.1) 

where j represents the polarization state either horizontal ( ^ ) or vertical ) 
polarization; K is the Beltzman's constant; V is the propagation frequency; c is 
the velocity of light, and J , S , 4* are the coordinates that specify the point 
of interest. 

In what follows we shall assume a uniform horizontally stratified atmosphere 
and choose a coordinate system with the y axis normal to the earth surface and the 
origin at the earth-atmosphere boundary directly below the point in the atmosphere 
under consideration (Figure 4.1). By expressing the apparent temperature at point 
■} and in the direction 0 and 4* as 4»), and assuming (1) the 

mechanism of raindrop scattering to be a random function uniformly distributed 
from 0 to 4 7r steradlan for both polarizations, and (2) the radiation process to be 
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Figure 4,1 The coordinate system in the atmosphere for expressing the 

incident and scattering of the polarized apparent temperature. 

dominated by thermal radiation, the radiative transfer equation may be written in 
the form (Stogryn, 1964) 

cos 0 ^ + I oC(}) -t- P^cCjj] (}■. ( J-3 

+ J y C <^°s(S>) Ta.^C}'^ ©s^ 4>s) S.n d % ci 

I an V ^^ 2) 

where TxjU , 4>) is the apparent temperature with ^ polarization at an altitude 

of 1 and in 'the direction 6 and <f>. . The first term on the righthand side 

of (4.2), namely o^Q^TairtT), 'S source function due to thermal radiation. ocQj 
is the total absorption coefficient of the atmosphere, T.ir C 3 ) is the thermometric 
temperature profile of the atmosphere. The second term on the righthand side of (4.2) 
is the source function due to particle scattering (for instance, ram drops) of the 
incident radiation from. all directions into the direction of interest. fscQ ) is the 

scattering coefficient of precipitation per unit volume at an altitude ^ ; y(cos®) 

is the scattering phase function of precipitation per unit volume with 


CJOS <S) cos 0 cos 0s ® C^S C4-4 >s) 
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(4.3) 



€)y 4^ = observing angles 
= scatferlng angles. 

A defalled analysis of ) and y (;cos <fi)J will be presenfed In fhe following 

section. * 

4.2.2 Apparent Temperature Theory with ^5^ Neglected 

Under clear sky, cloudy sky and light to moderate rain conditions the effects 
of scattering by the atmospheric constituents may be neglected for observation 
frequencies less than 20 GHz (Porter, 1970). Under these conditions^ the source 
function in (4.2) consists only of the thermal radiation term and (4.2) reduces to 
the form: 

cose J""} — -f C}, e, (.3-3 (4.4) 

For uniformly straitified atmosphere bounded below by an isotropic surface, there is no 
azimuthal variation. For this case (4.4) becomes 

cos ©-^U.e3 -i- e) TaiK-tJ) ^(^3 (4.5) 

Equation (4.5) is a first order differential equation. It can be solved with 
appropriate boundary conditions. The solution for a downward looking radiometer 
situated at an altitude ^ above the ground has been found to be (Sfogryn, 1967) 

T.jU.e) = L cj.e) [e.ce^T^ -t- T.-te)] 4- Ta+m(-3',e3 ( 4 ^^) 

j. =. ■ft. oi-* XT 

where apparent temperature at an altitude J , and nadir angle 0 

= transmittance ^ 

= exp sec ©j du. j (4.7a) 

oCiSj ) = total absorption coefficient of the atmosphere at an altitude ^ . 
6^C9) = surface emissivity (see explanation in Chapter III) 

= surface thermometric temperature 
atmospheric temperature 
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“ sec 


6 ] exp[-spc©j o(iLL:)au] d.^' 

Jo -'y 


(4.7b) 


'T^rU) = thermome^ric femperofure profile of the dir 
X,(») = reflected sky temperature 


= ^ f[\ Ce C6. ) s,. 6, d 6s d t. 


(4.7c) 


V; C6, ) = polarized differential scattering coefficients of the surface . 

d 

~^^(^) = sky temperature 

r V 

= sec Qs. Ta;^C >0 o<(J'; exp[-scc ©jJ << Lu) daj ot (4^7d) 


As shown in (4,6) and (4.7) the apparent temperature contains the attenuated 
surfoce brightness temperature, €^(9)T^ , os well as the direct and reflected 

temperature contributions from the atmosphere. To compute the apparent temperature 
Taj (>, &) f we use fhe surface roughness model developed in previous chapters 

0 Q 

and assume sea water .temperature to be 290 K. To account for the temperature 
contribution from the atmosphere the ARDC standard atmosphere which describes the 
temperature, the pressure and the water vapor profiles and the absorption coefficients 
of various constituents due to different weather conditions are needed. The weather 
model proposed by meteorologists that describes the spatial distribution of liquid 
water content and precipitation rates needed in the apparent temperature computation 
will be discussed separately in the following section. 

For an ARDC standard atmosphere, the temperature, pressure and water vapor 
profiles with altitude are given by: 


T„vt3-) = Tj - 6.S ^ -for ^ U 

= Tj - 7/. .5* •!< J >- l«n 

P tj ) = Pj exp F- O- 5 S’ ] 
f L> i r- S/2-2] 


(4.8a) 

(4.8b) 

(4.8c) 
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where. 




Tg = sea surface temperafure In • 

Pg = sea surface pressure In mm of Hg. 

2 

A = sea surface wafer vapor denslfy In grr/ m . 

^ = alfifude In km. 

For mulfiple consflfuenf atmosphere and by linear superposition, the total absorption 
coefficient c<(J) In (4.5) and (4,7) can be written in the form 

o^c£a(5-) + (4.9) 


where 

°tsU)= c lear sky absorption coefficient 
= cloud absorption coefficient 
= rain absorption coefficient 

If there is no cloud or rain, the total absorption coefficient is due to clear 

sky absorption only. The same kind of argument can be applied to the other weather 
conditions . 

A) Clear Sky Absorption Coefficient 

The two major contributors to the clear sky microwave absorption spectrum 
are oxygen and water vapor, thus, the clear sky absorption coefficient is 




"•f'Tkn 


(4.10) 


A.1) Molecular Oxygen 

The oxygen molecule absorbs microwave energy through the interaction of its 
magnetic dipole moment with the electromagnetic field. This phenomena was used 
by Van VIeck (1947a) to develop expressions for the absorption spectra of molecular 
oxygen in the microwave region. Oxygen has a series of absorption lines near 60 
GHz (often referred to as the oxygen complex) and a single line at 118.75 GHz. 

Using 60 GHz as a central line. Van VIeck derived an approximate expression for 
the absorption by the oxygen complex valid for frequencies outside the 50-70 GHz 
region. The dependence of the oxygen absorption on temperature and pressure was 
further investigated by many investigators (Hill and Gordy, 1954; Artman and Gordon, 
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1954; ZImmerer and Mizushima, 1961; and Meeks and Lilley, 1963), The more 
recent’ resulfs for calculating tlie microwave absorption by oxygen (outside the 
50“70 GHz region) may be summarized by the following expressions; 

(4.11) 


where. 


‘fco 




1 








( 6 o- ir^ ( 6 o + 7tJ2-|. ^^2 A 


7.]w- 


1 2a) 


°<n8 (?)= y 

^ - frequency in GHz 

83.50 

\^»i= 1. 20 


A 7^ 




(ug.TS -V vy -f 


] ,4. 


1 2b) 


V = 02 1 Pty) 

^ [TavC^)]^ 


= 1. 5 


f mi-\ /_113_V 
^ T6oM TaVii-)/ 




(4.13) 


P(l) = afmospheric pressure in mm Hg at altitude T". 

T^lr('3) = atmospheric temperature in °K at altitude . 

A .2) Water Vapor 

Though water vapor has only two absorption lines in the microwave part of the 
electromagnetic spectrum, at 22.23 GHz and 183.3 GHz, the cumulative effect 
of the large number of submillimeter water vapor line can be significant.. Using 
a low frequency approximation for the line shape factor. Van VIeck (1947b) derived 
an approximate expression for the residual effect of all the major water vapor 
lines whose resonances are higher than 30 GHz (thereby including the 183.3 GHz 
line). These expressions were later modified by other investigators (Bartett and 
Chung, 1962; Chung, 1962; and Staelin, 1966) to include the dependence on 
temperature, pressure and water vapor density using experimentally determined 
values for the line with parameter, Thus water vapor absorption can be 

expressed as: 

(4.13) 
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4 - ^ - 


( 22 . 235 -+ v-y -+ 


[t..v(^)J"^ 


"'P / Icr 


(4.14b) 


kiz. = 43.2 


-3 


^i-es = 2.55 X 10 

y„(fj= mif(?)( I -ho.li ^'*Vgy‘'^ )lT..v(H]'‘ ’^p[- 6 .VT...qj]( 4 . 14 c) 

Ay = 2.62(1 +0.001 ) F12 )/t 6 o CTi.+«)/s '8 ( 4 . 14 d) 


0 

= wafer vapor densify in gnr/m af alfifude }'• 

B) The Absorpfion Coefficient of Cloud ond Rain 

Ryde and Ryde (1946) were the first to calculate the attenuation of clouds 
and rain based on Mie*s theory (190Q, and their result has been the basis for much 
of the subsequent works. By using Laws and Parsons( 1943) drop“size distribution of 
various precipitation rates, they calculated the theoretical attenuation coefficient 
of rain directly in terms of eight different precipitation rates with eight wavelengths 
from 0o3 to 10 cm at 18^C. For temperature different from 18^C, the correction 
factor is tabulated from 0 to 40^C at lO^C interval. 

For non“precipitating clouds with drop radii much less than the wavelength, 
the familiar Rayleigh scattering equation was used by Goldstein (1951) to calculate 
cloud attenuation. 

Gunn and East (1954) not only reviewed Ryde's development since 1946, 
but added the backscattering cross section calculation for clouds and rain which 
has been widely quoted In radar meteorology. The dielectric constant of ice is 
very much different from water. Consequently they calculated the attenuation of 
ice and water clouds separately, and the temperature effect was also Included. 

Medhurst (1965), using Ryde*s theoretical approach, re-calculated the 
attenuation for rain with five precipitation rates and sixteen wavelengths from 0.3 
to 15 cm at 20^C. He also collected rain attenuation measurements from 1946 to 
1962, and compared with his theoretical result. He used a set of theoretical 
maximum and minimum attenuations for rain to compare with measurements and 
concluded that the agreement was not entirely satisfactory . There was a tendency 
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for measured a^^enua^^on to exceed the maximum possible levels predicted by the 
theory. But recent measurements by Browri and Haroules (1969) at 8 GHz to 15 GHz 
indicate that there is no tendency for measured attenuations to lie above the 
theoretical maximum values. Hence theoretical calculations seem to provide a 
reasonable basis for the prediction of attenuation by rain. 

Benoit (1968) interpolated Gunn and East’s discrete data with respect to 
frequency and temperature and have come out with one simple empirical equation 
very convenient for computation purposes, particularly when the temperature varies 
continuously with respect to height in the atmosphere. 

Following Benoit’s result, the absorption coefficient for clouds Is expressed 
in the form 


M 4^Ap^a.; (4.15) 

3 

— liquid wafer content- of clouds in grr/m 

V = propagating frequency in GHz 

t> = frequency index 

H-l. 

b„, = 1.95 for water cloud 

■ 1 ' 

bi, = 1.006 for i.ce cloud 

a = temperature coefficient 

- 6.-« Vfc C' (Ta.r(J) - 273jiJ ' ’ 

CLC, = - 8 . 26 1 1^1. - 1-767 X \o'^ ( - ZViJ)- 4,374- X Ta!rl}J - 

Figure 4.2 shows Gunn and East’s original discrete data and the straight 
lines were drawn by Benoit using the Interpolating equation (4.15). 

Presently available absorption coefficients for rain are all expressed in terms 
of attenuation in db/km. If the scattering effect is Ignored, the attenuation is 
mainly due to absorption and we can use attenuation coefficient to calculate the 
apparent temperature. Here we use Gunn and East’s empirical expression: 


(>)==. ;.a 


(4.16) 


to relate the attenuation to the precipitation rate. Is the attenuation in 

• . ! ' • ■ . * ^ ^ 
db/km and p is the precipitation rate in mn/hr. The parameters a and b are 
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frequency dependenf consfanfs. Referring fo Ippoilfo (1970), ihe a and b 


parameters are obtained 

with Laws and Parson's distribution at 

three frequencies; 

a = 0.008 

b = 1.32 

at 8.9 GHz 

a = 0.0125 

b = 1.25 

at 11.1 GHz 

a = 0.026 

b= 1.18 

at 13.9 GHz 


Ippolifo has shown fhaf the attenuation calculated by this formula is between 
Medhurst's (1965) theoretical maximum and minimum attenuation. 



Figure 4.2 Attenuation Coefficient for Clouds. 
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4.2.3 Apparent Temperofure Theory with included 

Under heavy precipitation a large percentage of precipitating particles 
exhibit diameters comparable to wavelengths at frequencies in the range 10 to 15 
GHz. As a consequence scattering cannot be neglected in considering the transfer 
of radiation through heavy precipitation. The physical effects are clear. The 
radiation from the surface not only encounters absorption but also scatters in all 
directions on its upward course. The addition of "collision" on its upward course, 
as a consequence, attenuates the surface radiation and reflections more effectively. 
Similar physical process describes the propagation of the inherent radiation of 
the atmosphere. 

To compute the apparent temperature 4*), the radiative transfer 

equation presented in the previous section should be employed and is rewritten 
below 

cos e j + fscC^^] CJ- ) 

2TTrTt 

J / C ©S d 9s d 

O ’^9 

^ or (4.2) 




t 


A close form solufion of (4.2) such os fhe soluMon of the first order differential 
equation presented in the previous section is unwiedly complicated and impractical. 

A finite difference method that simulates (4.2) as closely as possible at the net 
points is employed in this study. 

Lettjng /U. ^ cos ^ ^ == 

T — J (where t is called the optical depth of the atmosphere). 

We may express (4.2) in the form 

^ + TajC-r. a. 4 .) = Cl - 'CrCf “V- 


,1-n , I 


-V- -w ' 


(4.17) 


j VC C-OS ® ) Ta^ (. T Mj, *^4 ) 4*1 

5- -ft, or V 

with cos® ^ a/ I /U.4 co»sC4'-4‘i) (4.18) 

The direction cosine of M , M-s and the azimuthal angles 4> and is 
presented below (Figure 4.3). 
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Figure 4.3 The coordinate system using the direction cosine and the 

optical thickness for expressing the incident and scattering 
of the polarized apparent temperature. 


The scattering characteristics of precipitation represented by the function ^ , 

^ r ^ ^ f V C tPS 0) in (4.16) are also discussed below. 

Equation (4. 17) consists of three variables T , jlL , <p> . To further 

simplify it, the scattering phase function y ( Ccs , can be represented in a 
finite series of Legendre Polynomials (Shifrin, 1969). 

yccos®3= ^ v; p cc<,s®3 (4.i9) 


where Is ^he Legendre polynomial of order i . Addiflonal simplificafion will 

resulh when ^he following laenfil'y is employed. 

Ceos ®) = Pj -^/l J I -Ai CosC<t> - 4 .^ J 


= Pi (>* ) P^C/IO + 2 i P/ 0^3 CO s »iC4> - 4>5 3 

(4.20) 

where 

( ) = associafed Legendre polynomial of order ^ degree m . , 

We fhus have . , 


= y (m, 4>; , 


_ 

^ f\ Cos C - <f>s 3 
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(4.21) 



where 






(4.22) 






1 \ (^ TV\. = O 

O otUer va 3 ( s e 




EquaMon (4,21) suggests fhat we expand ^ , 4^ ) !n the form 

d 


N 


~ L t”<; ^ At) c 


>»i »o J 


• s y>\ 4* 


(4.23) 


When equations (4.21) and (4.23) are substituted into equation (4.17), the following 
independent equations result: 


^3 -T- JA C/t^^oT. (T,>tO 

— ( 0 


+ Cl - U> tr) ) Ta.v (T) 

C >^ = O, I, - ■ ^ 1sJ j , 


dMs. 
(4-.24) 


The boundary conditions for d rough surface bounded below and the atmosphere 
extended to infinity may be approximated in the form 


T 


. t r I - c-7.^ . 

C >t ) = o 


M. O 


^ O 


(4.25) 


where is the total optical depth of the precipitating atmosphere, e.{ /l) is the 

0 

rough surface emissivlty and Tj is the surface temperature. 

Let us assume (1) the rough surface is isotropic, and since (2) there are no 
sources above the atmosphere, and (3) the thermal radiation in the medium is also 
isotropic in nature. Then the solution of equation (4.24) consists of one term (m = 0) 
only and it becomes a function of T , /x , and independent of 4^ . 
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A) Mefhod of Finite Difference for Solving Radiative TransFer Equation 

In order fo obtain the function Tj(T, y(X ) numerically, the approximate 
problem is formulated by confining ourselves to (N = ) 5 terms in the expressions 
shown in eqatuions (4,21) and (4.23) and requiring that equation (4.24) and boundary 
condition (4.25) be fulfilled at the points /*,, ••• » interval (0,1) 

and n. , M. /•••i M , in the interval (0,-1). The integrals are here replaced by 
quadrature summations. 

A system of ordinary differential equations of the first order with N boundary 
conditions at t =0, and another N conditions at T = is obtained for the 2N 
unknown functions ^ = 1/2, . . .,N,~1,“2, . . .,-N) of this 

approximate problem. 

We introduce a difference net in t in such a way that the boundaries of the 
different substances in the points of discontinuity of the coefficients of equation 
(4.24) are points of the net. A finite difference system is obtained by integration 
over each step of the net (t^ , of the system of differential equations for the 
functions (‘^) on the assumption that these functions can be interpolated linearly 
over the step. For N = 5, the system of equations with the subscript j removed, is 
represented as follows: ' 


C' -^f^i)) T«,v CT^) -f 

, Tr S p 

2 )(Tsi,h ) -V (xt,7/Uc ) -h 




(4.26) 


^ = C ' - (r^) -v 
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)-^i(Tr + X| ) - c 

a. ^ ^ A, . Fa , .. „ \ \ . . 


2 ^ [A(a, •+t;^)' ) -+• AC/'v'^O^ V-+- 3] 


(4.27) 


^ T,:f ) ^ c ) = C ' - - cr.) ) T... Ct, ) H- 

a c [^ AO ( V-v T,-;, 3 -vA -Me )CV-v 1 ;:, 3 ] 


Note that 


i(M)dM *=_J [ .fOMj + f (-A)J £<^ 


t*' ? 


«. [ -fCAp i- 


(4.28) 


0| = 0.nS46, 02 = 0.23931, 03 = 0.28444, a^ = 0.23931, = 0.11846 

= 0.04691, ^2 ~ 0.23077, = 0.50000, = 0.76923, = 0.93509 

where d/ s are the weighting factors and A/^ are elements of the domain of 
T when a five point Gaussian quadrature integration technique is employed. 
Equivalently this can be written in matrix form as: 


®« V ■+ ( T,:, T,' ) = F, 

CAT,* -»- tA3 + A^T^" = F, 


(4.29) 
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(4.30) 


, 

C, = [ ccj ] 


. — ( I \ r K . , V 


^ ^ ^ X \ 'i 'TT OJ ( tQ * , , - 

.V 2 ^ 



Tr 

2 





^ = j- 

0+Vier\jD» se 


The boundary conditions of equation (4.25) imply: (a) on the rough surface 
T(r„) = ^ { Ci] - cf.-^ } . T'C'fo; 

where 



I. 6j IM,) . 


(b) with no incident radiation above the layer L 

T"Co) = [o] 

(c) an invariance principle is based on the general reflection ancJ emission properties 

of the layer , 




+ 'rr, 


(4.31) 


with 

^ |- -1 "t“ -f* 

cr, = [1] - ; -TT, e. 

where [ I ] is the unit vector, is the emission vector of the rough surface, and 
or^ , are the reflection and emission vector of the layer ( T^, r^). When 

equation (4.31) is substituted in (4.28) and solved for , and , we get 


[ C, - ce.,9^ + 13^,1 

<' [ C, - i- + /\0~ ' C f )] 

-t- c, )Cc,S.^A^)-'c, (4.32) 
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(4.33) 


If is also known from the fransmission’ properfies of fhe layer that 
where 

A 

^ = a^■^enua^^on coefficFenf of the layer 
= emission coefficient- of the layer 

From equations (4.33) and (4,28), the following expressions result 

^ - C AO"' 

S, = C 5. + A, r' [ F.:, - C, ( Tt,., T, )] <'*•“> 

The matrix functional equations shown in equations (4.28) to (4.34) which represent 
the original finite difference equation are applied to solve the integral differential 
equation with two boundary conditions given in equation (4.25), 

A A 

Compufafipnally, fhe coefficients ^ ^ are first found for 
-? = 0, 1, using o; = [ 1 ] ~ / and the recursive relations 

given in equations(4.32) and (4.34). These coefficients do not depend on 
With these coefficients calculated at every net point, equations (4.31) and (4.33) 
therefore enable us to find the solutions Tj^and for all values of ^ . To 
increase the accuracy of the numerical computation, the number of nets for and 
can be increased. Under the assumption that the emission induced by scattering 
contains equal probability for both polarizations, the system of difference equations 
presented in equations (4.28) through (4.34) Is applicable for either horizontal or 
vertical polarization; the only difference is the surface reflection and emission 
characteristics, namely o; = [ 1 ]~ ir^= Tj where j = -A or 2/ . 

B) Scattering Characteristics in Precipitation 

The microwave scattering characteristics in precipitation consists of the 
extinction coefficient, the scattering coefficient and the scattering phase function. 
Although the extinction coefficient has been computed for various precipitation 
rotes in radar applications (Ryde and Ryde, 1946; Gunn and East, 1954; Medhurst, 
1965), the scattering coefficient and phase function are still lacking. In what 
follows these coefficients will be derived first for a single drop and then the result 
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will be extended to a unit volume of water drops using on appropriate distribution 
function . 


(i) Single Drop Characteristics 

It will be assumed that the precipitation particles are spherical in shape so 
that Mie’s theory of scattering by a dielectric sphere is applicable. A detailed 
analysis of Mie's theory can be found elsewhere (Mie, 1908; Stratton, 1941; Van De 
Hurst, 1957). Here, only a brief outline is given. 

Suppose a plane wave E = ^ C ^ with unit strength, time 

dependence propagation constant , and polarization in the J direction 

is incident in the j direction (Figure 4,4), then the scattered electric field, 

due to a spherical particle situated at the 
origin, at large distance R is given in the 

A. 

spherical coordinates by 





E == 


= 


- e 


A (.9) Cos 






^ie) S(Vt 4>. 


(4.35a) 


(4.35b) 


Figure 4,4 The scafl'ered elecfric field due fo a spherical particle. 


where 


A (.6) 
3(6) 


^ + r - s FnCcosS j cl pjcccsff ) 

^ n^n-fO I S'le aie 

H « I ^ 

(Zk + f y r ^ ^ ^ pK ( <9 J ^ I s Pn ^ J 

o I st■>^ e J 

>c — I 


(4.36a) 


(4.36b) 


P. 


I 


C Cos G) 


associate Legendre polynomials 
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(4.37a) 


S __ S., ((•) . g.. Cm (>J 

^ <p^((> ) s^'c ~ ^ 4>/ (p) S^( 

g-.fr-’ ( p)^^('^p) 

S^C-yr,^)- i>^((>) (4,37b) 




-M 

^M. ( ^ ) = ^yiC^y 
C (X ) = C - ( )' 

(x) = 7^* order Bessel funcMon of fhe first- kind 
^ = normalized radius 

= 2T\a./x 
A = wavelength 

a. = particle radius (in the same unit as A ) 

= square root of the complex dielectric constant of the particle 
= (6'-^€")'"^ 

= f ( 1 - j;*^) 


The prime symbol in (4.37) indicates differentiation of the function with 
respect to the variable shown in the argument. 

The power scattered by the particle into a unit steradian in the direction 
( 0, t) is 


CC®, ± [|Aco3|^ cos^ 4> s-,w*q.] 


(4.38) 


For a randomly polarized beam, the scattered power may be obtained by 
averaging c ( © , ) over 'P , namely. 

US) ^ 

The scattering cross section )(( ) defined as the radio of total power 

scattered by the particle In all directions to the incident power can be computed 
by using (4.39) as 

f>) = 2TT ) CL9) SvM. 0 d 0 (4.40) 

■ o 
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The scaM'ering coefficien^ f^sc( ^ ^ furn, defined as fhe rafio of 

the scattering cross section to the geometric cross section of the particle, can be 
computed by using (4.40) and (4.36), It Is expressed in the form 


Similarly, the extinction coefficient is expressed as (Van de Hurst, 1957). 

I 

f )=, O Re C ^ 


(4.41) 


(4.42) 


The scatfering phase funcfion, 'f ( 9 ) , is defined as the ratio of the 

scattered power directed into a unit solid angle in the direction 0 to the scattering 
cross section . It is expressed in the form 


ue) _ I AIPJl^ -t- iBlejp 

2Tr<x^ f ) (4.43) 


From (4.43) it is obvious that • - - 

j 4 ce) si’vu e d© = T , 

Equation (4.41) through (4.43) are the formulas from which the scattering 
characteristics of .a single drop can be calculated. 

(ii) Unit Volume Characteristics ^ . 

Under the assumption that many particles of different sizes are randomly 
distribufed in space withlii a unit volume of precipitation, the macroscopic scattering 
characteristic can be computed by integrating the single drop characteristics with 
a known drop size distribution functlon over the entire size distribution range. The 
scattering characteristics of precipitation for a unit volume may thus be expressed in 
the form 

t • ' ■ w' ' ‘ ' 

= TT J Cl* tr ^ f ) N (a) 4 

» • . i . . . ■ 

tfO 

^sc TT J a" C f ) M (a)' d. a. 


(4.44) 

(4.45) 
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yc®) = 



H-'sc C -f N(o) ^ a. 


(4.46) 


where . 

Pex = extincflon coefficient per unit volume of precipitation in neper/cm 
^sc - scattering coefficient per unit volume of precipitation in neper/cm 
= scattering phase function per unit volume of precipitation 
^ = radius of raindrop in cm 

= drop size probability density function in the precipitation 
For the case of a Marshall-Palmer distribution, the probability density function is 
expressed in the form 

^l (a) = exp C- C, a ■) (4.47) 


where No = S x lo ^ ^ 

- D. ^ 1 ‘ 

C. - r om-' 

■p = precipifaHon rafe in millimefer per hour 

(Hi) Computed Scattering Characteristics of Precipitation 

From the above derived equations, the extinction coefficient, the scattering 
coefficient, and the albedo of scattering can be computed and are given in Table 4.1 
with different frequencies, precipitation rates,, and raindrop temperatures i The 
scattering phase function versus angles are plotted in Figure 4.5 through 4.8. 

From Table 4.2, it is shown that extinction and scattering coefficients 
increase with increasing precipitation rate and, frequency . On the other hand a 
decreasing raindrop temperature will increase the albedo of scattering. For pre- 
cipitation rate less than 10 mn/hr and microwave frequency below 9 GHz, the 
scattering effect is negligible since the contribution of the scattering term is less 
than five percent of the total energy. 

From Figure 4.5 through 4.8, it is shown that an increasing precipitation 
rate at 13.9 and 8.9 GHz will increase the backscattering power. As mentioned 
before, this scattering phase function can be approximated by a finite series of 
Legendre polynomials 


•f 
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TABLE 4. 2 

THE EXTINCTION COEFFICIENT, THE SCATTERING COEFFICIENT, AND THE ALBEDO OF PRECIPITATION 
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4.3 COMPUTED APPARENT TEMPERATURES OF CLOUDS AND RAIN 


4.3.1 Infroduction 

The fheoreMcal invesMgafion given in The previous fwo secMons described 
in general fhe equations to be used for various weather conditions. To compute the 
downward- looking apparent temperature for a ^s^excluded weather condition, it is 
necessary to assume some values for the parameters in (4.6) and (4.7). As can be 
seen in (4.6) and (4.7) both sea surface brightness temperature and atmospheric 
contribution play an important role to the total apparent temperature. Here we 
employ the results of Chapter III and a sea water temperature of 290^K to calculate 
the rough sea surface emissivity at 13.9 GHz. Since sea surface roughness charac- 
teristics can be related to wind speed over the sea as mentioned in Chapter II, in 
what follows we also assume a wind speed variation of 5 to 15 iVsec to represent 
different sea surface roughness. 

To compute the appaent temperature due to atmospheric contribution it is 
necessary to assume the atmospheric temperature profile, and the absorption 
coefficient for a specific weather condition. In this study unless otherwise stated we 
apply ARDC standard atmosphere for describing atmospheric temperature profile, 

. To compute the clear sky absorption coefficient we also employ ARDC 
standard atmosphere in which the temperature, pressure, and water vapor profile of 
atmosphere are given by (4.8). The magnitudes of temperature, pressure, and water 

vapor density at. sea level are yet to be specified. In the following, an air temperature 

o 3 

of 290°K, a pressure of 760 mmHg and a water vapor density of 7.5 g/rn at sea level 

ore assumed. 

The presence of clouds and light to moderate rain can be expressed as 
additional absorption to the already present clear sky absorption and the ARDC 
standard atmosphere is still applicable. But the spatial distribution of the liquid 
water content for clouds and the precipitation rate for rain are required to compute 
the cloud and rain effects. Some realistic clouds and rain models proposed by 
meteorologists which include the spatial distribution of liquid water content for 
clouds and precipitation rates for rain are presented below. Subsequent computations 
of the apparent temperature using these models assumes an observational height of 
7 km. 
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APPARENT TEMPERATURE 
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FIG. 4.9 The Computed Apparent Temperature for Observation 
Throuyh Stratus Clouds 

7l 


4.3.2 Clouds 


According to meteorologist classification (Donn, 1964) there are three basic 
cloud forms^ namely cirrus, cumulus, and stratus. All cloud formations consist of 
these standard forms, a combination of these forms, or a modification of these forms. 

Since cirrus clouds are typically thin layers of ice, their absorption (or 
emission) are relatively insignificant compared with the other cloud types at micro- 
wave frequencies. Hence only stratus and cumulus cloud formations are considered. 
Usually the liquid water content of the cloud varies with height and attains its 
maximum near the cloud top. However, for convenience of computation, the average 
value over the total cloud thickness is employed in the calculations to follow. 

(i) Stratus 

Stratus clouds are typified by a thin layer of low water content clouds. The 
apparent temperature versus nadir angle calculation of these cloud conditions are 
based on Neiburger*s models (1949) given in Table 4.3. 

TABLE 4.3 

/ 

Neiburger*s Stratus Models 


Cases 

Altitude Extent 

Water Content 


(m) 

(grn/m3) 

Case 1 

30-580 

.35 

Case 2 

152-520 

.25 

Case 3 

213-490 

.20 


The results for horizontal polarization at 13.9 GHz and 5 rr/sec wind speed 
are shown in Figure 4.9. 

It is evident that stratus do not contribute much to the rise in the apparent 
temperature compared to the clear sky case. 

(ii) Cumulus 

Cumulus clouds are prevalent over the ocean for all seasons. To typify their 
effect Levine's models (1965) for cumulus are employed (See Table 4.4 below). 
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FIG. 4.10 The Computed Apparent Temperature for Observation 
Through Cumulus Clouds 
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APPARENT TEMPERATURE 
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FIG. 4.11 The Computed Apparent Temperature for Observation 
Through cumulus Clouds 
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TABLE 4.4 


Levine's Cumulus Models 


Case 

Altitude Extent 

Water Content 
(grn/ m'^) 


(m) 

Case 1 

457-1068 

0.5 

Case 2 

. 457-2590 

1.0 

Case 3 

457-3810 

1.25 


The computed apparent temperature versus nadir angle by using (4.6) and 
(4*7) for these three cases are shown in Figures 4.10 and 4.11 for horizontal and 
vertical polarization, respectively. Results are also shown for different surface 
roughness conditions as well as cloud cases. The result indicates that (1) for the 
some surface condition the thicker and the higher the water content, the higher the 
temperature rise with respect to clear sky; (2) at larger nadir angles say around 50 
degrees horizontal polarization is more sensitive to the temperature rise as compared 
to vertical polarization; (3) the rise of apparent temperature ‘due to the increase of 
surface roughness, namely increasing wind speed, can be detected in all cumulus 
cases and (4) comparison of the horizontal apparent temperature curves between clear 
sky and cumulus cloud (Case 1) (Figure 4.10) a qualitative differentiation between 
surface effect and cloud absorption effect can be made. The cloud absorption 
temperature increases with increasing nadir angles, while the surface effect causes 
the same amount of temperature increase over all nadir angles. 

(iii)’ Overcast 

Overcast may be considered as a combination of stratus and cumulus clouds. 
Frequently an overcast condition can be associated with the warm sector preceding - 
a cold front. An overcast may be classified according to its depth and water content. 

To study these overcast conditions Porter's models (1969) are employed. A 
description of these overcast situations is shown in Table 4.5. 
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FIG. 4.12 The Computed Apparent Temperature for 
Overcast Conditions 
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FIG. 4.13 The Computed Apparent Temperature for 
Overcast Conditions 
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TABLE 4,5 


Porfer*s Overcasf Models 


Classification 

Altitude Extent (m) 

Water Content 


(gnr/m^) 

Light (sun visible) 

300-650 

0.33 

Medium (1 ight sky) 

400-900 

0.67 

Heavy 

500-3200 



The apparent femperafure characl-erisl-ic for these overcast conditions are 
shown in Figures 4.12 and 4.13 for horizontal and vertical polarization with sea 
surface wind speed of 5 rr/sec and microwave frequency of 13.9 GHz. In the figures, 
the horizontal polarization exhibits an increasing sensitivity with nadir angles under 
the overcast condition. A 70^K separation between the heavy overcast and the 
clear sky condition can be observed at nadir angles in the vicinity of 70 degrees. 

The response for the vertical polarization on the other hand shows maximum sensitivity 
to clouds only at nadir. 

4.3.3 Light to Moderate Rains 

The extent of precipitation over the world at any time is small in comparison 
to the extent of the non^precipitating regions. The frequency of rain over oceans is 
known to be even lower than that over land. As a consequence the operation of a 
combined radiometer-scatterometer will not be impeded often by rain. It is important, 
however, to estimate the effects of rain on the radiometric observation to attempt to 
discriminate between precipitating and non-precipitating regions. 

The characteristics of precipitation have largely been described over land. 

No information regarding vertical rainfall distributions has been found for rain over 
oceans. It Is assumed that the rain characteristics over land and sea are similar. In 
this respect Valley's rain models (1965) are used in computing the apparent microwave 
temperature. He described several cases, only one case, however. Is treated here. 
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FIG. 4.14 Precipitation and Water Content Distribution 
for Extensive Showers 
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FIG. 4.15 The Computed Apparent Temperature for Observations 
Through Rain 




APPARENT TEMPERATURE 



This case is rcpresenfal ivc of summer rain in temperate latitudes and with a slight 
increase in temperature, of widespread tropical rains as well. 

This model typifies widespread uniform precipitation and not showery 
conditions. The precipitation and cloud characteristics are shown in Figure 4.14 
and Table 4.6 respectively. As a result of different updrafts, four different 
precipitation rates can be described for each case. 

TABLE 4.6 


Valley's Rain Model 


Updraft 

i 

Condition 

(rr/sec) 

Rain 

Parameters 

Cloud Parameters 

Altitude 
Extent 
• (m) 

Precipitation 
at z = 0 
(mrr/hr) 

Altitude 

Extent 

(m) 

Water 

Content 

(grr/m^) 

0.4 

0-3100 

10.3 

3100-7000 

0.30 


0-3200 

7.9 

3200-7000 

0.25 

mg^ 

0-3300 

5.2 

3300-7000 

0.15 

0.1 

3-3500 

2.8 

3500-7000 

0.10 


The computed apparent temperature characteristics are shown in Figures 4.15 
and 4.16, From Figure 4.15 it is evident that the sensitivity to rain is greater at 
higher angles for small rainfall rates than for larger rainfall rates. In Figure 4.16, 
an approximately linear correlation of apparent temperature rise with respect to 
precipitation rate are demonstrated at nadir. From this it may be possible to predict 
precipitation rate from the measured apparent temperature whenever the surface 
characteristic is known. It is also clear that the radiometric contributions by rain 
are notably larger than those of clouds. 

4.3.4 Heavy Rain 

Based on finite difference method described in earlier sections, computer 
programs are developed to compute the apparent temperature as observed at altitudes 
above heavy rain. The results for two uniformly distributed precipitation rates, 10 
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FIG. 4.17 The Computed Apparent Temperature for a 
Precipitation Rate of 10 mm/hr. 
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FIG. 4.18 The Computed Apparent Temperature for a 
Precipitation Rate of 10 mm/hr. 
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FIG. 4.19 The Computed Apparent Temperature for a 
Precipitation Rate of 30 mm/hr. 
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FIG. 4,20 The Cbniputed Appare^nt Temperature for a 
Precipitation Rate of 30 mm/hr. 
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wrr\/hr and 30 mnr/hr, are shown in Figure 4,17 through 4.20 for both horizontal 
and vertical polarizations. 

In the figures both scattering'included and“"excluded cases are given for 
comparison. For scattering-excluded case, we assume the albedo of’ scattering 
to be equal to zero in the radiative transfer equation given in (4.33), Since the 
source function induced by the scattering term Is not as effective as the thermal 
radiation source, a decrease of the computed temperature Is clearly evident when 
the scattering effect is included. At larger nadir angles (greater than 60^) the 
effect of scattering is seen to increase and causes a larger decrease in the observed 
apparent temperature. The effect of scattering is most significant as the precipitation 
rate increases. 


4.4 CORRELATION OF EXCESS APPARENT TEMPERATURE WITH ATTENUATION 
Introduction ’ 

The primary advantage of having a radscat system as opposed to a scatterometer 
system is that it allows correlation of the excess apparent temperature (i.e., the 
difference between the total apparent temperature and the surface brightness 
temperature) with the attenuation, thus providing a correction factor for the 
scatterometer data. This, in turn. Improves the wind speed estimation over the 
sea when cloud and rain intervene. 

Direct numerical computation employing (4.6) to plot the attenuation given 
by versus apparent temperature difference is possible but tedious. In what 

follows a simplified model is used to investigate the correlation between the excess 
temperature and atmospheric attenuation. 

The apparent temperature for a downward looking radiometer located at an 
altitude j above the sea surface with nadir angle 0 , as derived in section 4.2,2 
for non“Scattering atmosphere, has the form 

Ta 


+ T,^Ce)] 


(4.6) 


1 = 
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where J refer ^o fhe sfafe of polarization; L ( 0 ^ 3 ^ is the transmittance; 
is the surface brightness temperature; temperature emitted by the 

atmosphere between the surface and the height in the direction 9 
is the reflected sky temperature, Tr^ 19} is also expressible in the form 

Trj'®) = :fi;J V: (e, 9s, <PO'T=k)(9*) (4.7c) 


r.j (e ) 


w 


ith 


2-H 

G-^e)= \ J S''-9^d9sd4>s 

3 0 o 


^ is either or V" and H ^ 5 ^ ^s 3 is the surface differential scat- 

9 ^ 


where 

tering coefficient. 

Due to the' complexity of ©5 <)=»s) os derived in Chapter 3 and the 

integration needed to obtain the reflected sky temperature, the problem can only be 
evaluated numerically. However, in numerical computation the parameters that 
affect the apparent temperature cannot be easily visualized. To circumvent this 
difficulty a simplified model is employed to analyze the correlation of the excess 
temperature with attenuation. 

Let 








i 1 V®' , (4.48a) 




VT 




Y (6, 0s, %) 

Q 


(4.48b) 


== ' - ^(9) (4.48c) 

t 

where represents a weighted average sea surface reflectivity with respect 

to specific sky temperatures. Under many circumstances involving clear sky, stratu^s, 
cumulus, overcast and light rain at frequencies below 15 GHz, and nadir angle 
less than 70 degrees, it is found that insensitive to Tst^ W variation or 

atmospheric changes. It mainly depends on surface characteristics as can be 
expressed in terms of (4.48b) or (4.48c) without significant error. 
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Using (4,48) , we can express (4.6) in the form 


T,. I.G, n = L IS. 3 ) [ Tj ^ Xtjte ) ] + ,3 ) (4 . 


49 ) 


The significant advantage of (4.49) over (4.6) is that the reflected sky temperature 
'Ti.iCB) is now expressed as the product of the effective surface reflectivity fxtO) 
^C6-) the sky temperature, and that depends only on surface 

characteristics. 

To further simplify the problem, can also be written in the form 

(Peake, 1969) 


TsfejCe) ^ {\ \2 - 5o)( \ - e ) 


(4.50) 


where is the air temperature at the ground level; and is the total 

optical depth at nadir. The total optical depth is given by 


= J 


(4,51) 


where ) is the absorption coefficient of the atmosphere per unit length at 
altitude J . 

Similarly, we can also express ) >n the form 


— X Sec © 

- ^o) ( ^ - e ) 


wifh 


and 


X = 


UG, > ) = e 


T SGC 0 


(4.52) 

(4.53) 

(4.54) 


By employing (4.49) through (4.54), the' correlation of excess temperature 
with attenuation is presented below. 


4.4.2 Analysis 

To make the analysis simpler we assume the radiometer to be located high 
enough that the atmospheric temperature 3 ) is about the same as the 
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sky temperafure, i.e* 
where 


v/- •^roStc-Q\ 

T«-l»,3>= T„,W) = •’ (4,55) 

= T„ H 1 - L ] 


T„ = - 50 

, —To sec 9 . 

L = e 

with 3 so large that we may assume T = To wil-hout significant error. Using 
(4.55), we can express (4.49) in the form 


Taj = L [e^(.e)Tj -v T,., ( l - L )] -a-T^ ( i - Lj 


'3 


(4.56) 


When we define the excess apparent temperature, , as the difference 

between the total apparent temperature and the brightness temperature due to sea 
surface only, we get 

= c > - L)T„ [ 1 ^ 1- fj»)] + L 

Letting the attenuation factor. A, equal to the inverse of the transmittance, namely 

Sec 6 To 


ft = = e 


(4.58) 


K = See 9 -To 


We get 




(4.59) 
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ATTENUATION 











From (4.59) and (4,60) some remarks can be made: (1) For fixed Tj , and 

(and hence (0) ), |-he larger ^he oM-enuaHon A fhe larger fhe excess 

apparent- femperafure, , and fhe rafe of increase of fhe excess apparenf 

femperafure / decrease wifh fhe increase of fhe affenuafion A. (2) For 

a fixed affenuafion A and surface femperafure T*j , decreasing surface emissivify 
C6 ) will increase bofh fhe excess opparenf femperafure and fhe rafe of increase 
of fhe excess apparenf femperafure with respecf fo affenuafion. 

The frend of fhe excess apparenf femperafure versus affenuafion is demonsfrafed 
in Figure 4.21 . 

In Chapfer III if has been shown fhaf fhe surface emissivify decreases wifh 
increasing nadir angle. Thus fo obfain a beffer sensifivify of fhe excess femperafure 
versus affenuafion, if is suggesfed fhaf fhe radiomefer should look af larger nadir 
angles, say 50 degrees or larger, and operafe in horizonfal polarizafion mode. 
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CHAPTER 5 


SUGGESTIONS FOR FUTURE STUDIES 


The non-coherenf composite surface model is shown to predict scatterometric 
and radiometric responses of the sea with reasonable accuracy. The use of an 
isotropic surface model, however, has completely ignored the anisotropic character- 
istics of the sea surface. The cross wind response, for instance, cannot be distin- 
guished from that of the upwind. Also, the evaluation of the scattering integral by 
the stationary phase approximation results in a wind speed dependence of the large 
undulations only through the rms slope of the large structures and not the specific 
shape of the surface spectrum. It is suggested that future effort for sea scattering 
theory should accommodate the two dimensional (anisotropic) sea spectrum proposed 
by oceanographers. This approach may, however, be too involved for microwave 
emission theory due to the many lengthy numerical integrations required for a 
number of multiple integrals, but is certainly feasible for radar backscattering with 
the present computer capabilities. 

Theoretical computations of the microwave apparent temperature over the 
ocean have shown a definite temperature increase due to clouds and rain. However, 
the cloud and rain models presented are not extensive enough to include all possible 
conditions. The horizontally non-uniform cloud and rain models that will occur a 
lot of times should be considered in future studies, although it may further complicate 
the formulation of the problem. For heavy rain, the solution is based on the uncoupled 
radiative transfer equation. It may also be extended to a more complicated coupled 
integral differential equation as given in Stogryn's work. 

A simplified model is used to study the parameters that affects the apparent 
temperature under cloudy and rainy conditions. Good agreement between simplified 
model and detailed numerical computation for correlating excess temperature with 
attenuation suggests that the simplified model can also be employed for the sensitivity 
study of surface parameters such as temperature and salinity of sea water. Further 
correlations of the excess temperature for detecting a definite rain fall rate in a 
precipitating atmosphere may also be possible if the related parameters can be 
estimated properly. 
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CHAPTER 6 


CONCLUSIONS 


The invest igaHon in the previous chapters presents a complete microwave 
apparent temperature theory over the ocean. The model for computing the sea 
surface emission characteristics is the same as the model for deriving the expressions 
of differential scattering coefficient, since the surface emissivity depends on the 
differential scattering coefficient of the sea in an integral form given by Peake. On 
the other hand, the rise of the apparent temperature over the sea due to clouds or 
rain depends on its absorption (emission) characteristics, atmospheric temperature 
profile, and the specific meteorological model used. 

A bistatic two scale non~coherent scattering theory extended from Symyonov*s 
paper has been developed to yield the expressions for the differential scattering 
coefficients. The emission and the backscattering characteristics are then derived 
from the differential scattering coefficients in the standard way. The theory 
assumed Gaussian surface height distributions and Gaussian correlation functions for 
both scales of roughness. 

The emission and the scattering characteristics are shown to be dependent 

on the rms slope of the large undulations m, the standard deviation of the small 

irregularities cTj , and the correlation length of the small irregularities i . The 

wind dependence of the first two parameters is associated with m through slick sea 

measurements by Cox and Munk, and the through the high frequency sea spectrum. 

The parameter i is associated with the shape of the high frequency sea spectrum 

“4 

and can be reasonably chosen by fitting the sea spectrum BK to the assumed 
Gaussian spectrum. It is noted that the emission characteristics for horizontal 
polarization Is a sensitive measure of c?; , Thus, is established by fitting the 
emission characteristic to measured data for different wind speeds. The parameters 
chosen in this way are then used to compute the vertically polarized emission 
characteristic. Good agreement with measured data and better agreement than single 
surface model are demonstrated. 

The same set of surface parameters at each wind speed is then used to 
compute the backscatter characteristics. The results except for level are shown to 
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agree reasonably over all angles with NRL backscatter data under similar wind 
conditions. Comparison of these characteristics with a single parameter surface 
model demonstrated better results. 

These findings have proven that the validity of the assumed scattering 
model is better demonstrated when both the predicted backscatter and the emission 
characteristics are compared with measurements. They have further shown that the 
measured emission and scattering characteristics with the aid of a reasonable 
composite surface theory may aid the oceanographer in identifying the wind 
dependence of the sea spectrum. 

Numerical technique has been applied to solve the radiative transfer equation 
and to determine the cloud and rain effects on the observed apparent temperature 
when looking downward over the ocean. For large rainfall rates the effects of 
scattering are found to be significant and are incorporated into the radiative transfer 
equation . 

It is shown that stratus, cumulus, overcast, and rain all contribute significantly 
to the observed terriperature . Larger sensitivities to clouds and rain are observed for 
horizontally polarized apparent temperature at large nadir angles than for vertically 
polarized apparent temperature. 

Aside from roughness detections, the apparent temperature of a radiometer 
can be used to infer sea water temperature when surface roughness parameters is 
being measured by a scatter'ometer. 

The rise of apparent temperature has been found to correlate with the atmo- 
spheric attenuation. This result is useful for correcting scatterometric observations. 

It may also be used to differentiate between a precipitating and a non-precipitating 
atmosphere when the surface condition is known. 
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APPENDIX 

FORTRAN IV COMPUTER PROGRAMS 


A lisf of FORTRAN IV computer programs and a brief description of the input 
ond output format are presented in this appendix. The equations and formulas shown 
in the programs are based on the derivations given in the main text with some 
modification. The tabulated results and the figures of the main text are computed by 
applying these programs. 

A-1 A computer program for computing the sea water dielectric constant and the 
plane sea surface emissivity . 

In this program the input cards are punched in data”card format with different 
frequency FR(I), temperature TEMP(I) and salinity SALI(I), The output for the 
computed dielectric constant is tabulated with respect to temperature and salinity. 

A different frequency results in a different table. The output for the computed plane 
sea surface emissivity is tabulated with respect to nadir angle at a ten degree interval 
for both horizontal and vertical polarizations at one frequency, one temperature and 
one salinity. For multiple frequency, temperature or salinity input, the output of 
the computed emissivity is printed out sequentially with preassigned sequence until 
all values are printed out. 

A“2 A computer program for computing the surface emissivity due to small scale 
differential scattering coefficient . 

This program is written for computing the surface emissivity due to small scale 
roughness. From this program the contribution of the small scale roughness to the 
overall emissivity can be estimated. The inputs are punched In constant^card format. 
The first three constant^cards contain the information on surface roughness parameters, 
i.e. the rms slope of the large scale sea surface, and the correlation distance and 
the standard deviation of small scale roughness. On the second two constant-cards 
are the real and the imaginary parts of the dielectric constant. The output of the 
computed surface emissivity is printed out versus nadir ongle at a ten degree interval 
for both horizontal and vertical polarization. To compute a different surface 
condition either due to surface roughness or dielectric properties It is required to 
change the Input card. 
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A~3 A compufer program for compufing the modified Fresnel reflccHon coeffic?en^ , 

In fhls program fhe inputs are punched in constant-card format. The first fwo 
constant-cards contain the correlation distance and the standard deviation of small 
scale roughness. On the second two constcn t-cards are the real and the imaginary 
parts of the dielectric constant. The output of the computed results are printed out 
versus nadir angle at a fifteen degree interval for both horizontal and vertical 
polarizations. In this program Romberg integration is being employed. 

A"4 A computer program for computing the two-scale rough surface emissivity . 

The computed results of the last two programs and the differential scattering 
coefficient due to the large scale roughness have been employed in this program for 
computing the two-scale rough surface emissivity. Due to the complexity of multiple 
integration and the averaging procedure only one case is computed in one program 
time. The inputs are all punched in constant-card format. The first three cards 
contain the surface roughness parameters. The second two cards contain the real 
and the imaginary parts of the dielectric constant. The output of the computed surface 
emissivity is printed out versus nadir angle at a ten degree interval for both horizontal 
and vertical polarizations. 

A-5 A computer program for computing two scale backscattering cross^-section . 

In this program the input is punched in data-card and constant-card format. 

The data-card co'ntains the surface parameters i.e. the large and the small scale 
roughness. The fwo constant-cards contain the real and the imaginary part of the 
dielectric constant. The output of the computed backscattering cross-section is 
printed out versus nadir angle at a ten degree interval for both horizontal and vertical 
polarizations . 

A-6 A computer prdgr.am for computing the clear sky and cloudy sky apparent 
temperatures . 

In this program the Input data contains three parts. The first part is in read- 
card format that reads in the required surface emissivity and averaged surface 
reflectivity. The surface emissivity data are obtained from previous computer 
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program. The second porf of fhe input Is in data-card format. It specifies the 
standard atmospheric conditions such as the temperature, the water vapor and the 
pressure at the sea level. The lost part of the input data is in read^card format that 
reads in various cloud conditions, i.e. the cloud thickness and liquid "water contents. 
The output of the computed apparent temperature is printed out versus nadir angle at 
a ten degree Interval for both horizontal and vertical polarizations. Each cloud 
condition is printed out in one tabulated form. 

A“7 A computer program for computing the apparent temperature of light to moderate 
rain > 

The input and output format In this program is the same as the above clear and 
cloudy sky apparent temperature program. The only difference is in the third part 
of the input cards. In this program various rain cases are employed instead of cloud 
conditions. 

A"8 A computer program for computing the scattering characteristics of rain. 

In this program all the input data are punched in read-card format. The first 
card contains the number of integrations and the differential increment needed in 
the integrations to obtain good accuracy. The second card contains the number of 
frequencies, temperatures, precipitation rates and angles that are required in the 
computation. The last three cards contoin the values of the frequencies, the rain 
drop temperatures, and the precipitation rates. The output is headed by the title 
of the problem. For a specific case with one frequency, one temperature and one 
precipitation rate, the output of the computed results is printed out In two parts. 

The first part contains the extinction coefficient, the scattering coefficient, and 
the albedo of scattering. The second part contains the scattering phase function at 
nineteen different angles at a ten degree Interval. For different precipitation rates, 
temperatures and frequencies, the format of the printed output will repeat as above. 

A-9 A computer program for computing the apparent temperature of heavy rain. 

This program contaTns many matrix operations. The matrix Inversion sub- 
routine is obtained from computer subprogram library. Hence we only use a call 
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sfafemenf to call In this subrouHne. The inpul* data are punched in read-card forma^. 
The fiVsf card specifies fhe surface emissivity and reflecf ivify , The second card 
specifies fhe scaffering characferisfics of rain and fhe values are obfained from fhe 
above program oufpufs. 

The oufpuf of fhe compufed resulfs are prinfed ouf wifh respecf fo nadir angle 
and incremehf layer. In fhls program we assume fhirfy incremenf layers for fofal 
rain fhickne^s and five nadir angles. Hence, fhe prinfed oufpuf will confain one 
hundred and frffy nef poinfs. For differenf precipifafion rafes and surface condifions 
fhaf include fhe change of polarizaflon, if is necessary fo change fhe inpuf card and 
re-run fhe program fo obfain fhe compufed resulfs. 
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A“l. A computer program for computing the sea water dielectric constant 
and the plane sea surface emissivity . 


S IDFNt 7 57.3 WtJ SHIH TSENG 

t OPTION FORTRANtFRrNT/15/ ^ 

Ti FORTRAN Nr)Fr< 

DIMENSION FP(10,6,6) *FPP(10,6,6) »FR( 10) »TEMP(6) ,SALI (6) 
dimension FMH( 10 ,6 , In ) ,FMV <1 0 , 6 » 1 0 ) » T BH ( 1 0 ♦ 6 . 1 0 ) , T PV I 1 U » 6 » U> ) 

DATA (FR(I),I = 1,10)/1 .228,1.41,A.A55,8.36,«.91,9.3,13.V,15.6,17.‘. 
1.19,34/ 

DATA (TEMPI I). 1 = 1. 61/ 276. .2 8Y,. 284,. 288. ,292 . ,296. / 

DATA ( SAL I ( I ) , 1 = 1 ,6) /3 2. , 3 3. ,34. ,35. ,3 6. ,37. / 

DO 70 I =1 ,10 

RT^FTTm 

DO 70 J = l ,6 

TG^TFMPm ' : ^ ^ 

DO 70 < = 1 ,6 

S = SALI M< 1 ^ ^ 

CALL DTFLrO(FQ,TG,S,Fl ,E2 1 
hP'( I ,J,3i )=F1 
FPP(I,J.<)=F2 
/o continue 

DO 71 1=1,10 

W1TTTF~r5T7TT Firm 

72 FORMAT ( IHl ,//45X,5HFKE0= ,F8. 4, 4H GHZ,//) 

WRlittb, /3) TTFMPTTrrTJHTTFT 

73 FORMAT { 1 IX ,6( 5HTEMP=,F5. 1 ,2H K,4X)//1 

WRTTE(6",22 ) 

22 format ( •’X ,4H5AL I,6(4y,2HF1,6X,2HE2,2X)//) 

w(rnrr'6T74T rs’AunK) » (tPi itjti^) , tPF(-r7jT7:-r,j=r.6i';<=T.&) 

74 FORMAT! 1 3 ( 3XF5 . 2 ) / / ) 

rr Tom'mu'F — ^ 

PI = ATANM . )#4, 

P’lTD'i^ T 7T8TD ^ 

L = 4 

DU“H' T=TTI-n 

DO 8 J=1 ,6 

00" 8 ■< = ! ,i'0 

E1=EP( I , J,L ) 

^ F2 =EPP ( r; jTL ) ^ 

THFTA1 = FL0AT.( <-l )*PI10 
Ch S = cOs ( iHh'l A 1 ) 

SFN = 51N(THFTA1 ) ' 

ro.F7=rE'^'#rF.s 

<;i N? = SFN*‘^FN 



B=(Fl»Fl+r2*P?)**n.25 

G7SFrA O“n ,"5 » A I A N ( F2 /T i > 

F2=CFS*C0S(GAMA0) 

F = 1 , n / ( C 0 s 2 TT7(T5fTrrF7'+'Fti ID 

EMH ( I , J,K ) =4. »R»F2 »F 

GAmT2'^T7VN (T 27 1! H 

G=CFS<f ( F1*F1+F2^fE2 ) »»0. 5 







- KG=R*G*rO'; fG'AW7;^MAn ) 

R/4 = 1 ,0/ ( F1*R + G*G + 7 ,0*KG ) 

n^VCT *T»I^ ) =4. 'n?r,*R4 

TRH( r » J»IC)=TFMP(J)J»FMH( I ,J,K ) 

“T'W( I ♦T»l<:'T = TtMP(JI»FM7( I «J»K ) 

8 CONTINIJF^ 

snrmrrsTTTT 

21 FORMAT ( IHl »50HTHIS IS TO CALCULATE THE EMISSIVITY OF SEA SURFACEt) 

Do~ri r=i,io ^ 

WRITE(6,18) ( EMH( I ,5»<) ,K = 1 , 10 I » ( EMV( I »5,K ) »K = 1, 10) 

18 FORMAT ( 10I4XF 6.4)//) 

11 CONTINUE 

DO \2 J=1 *6 ^ ^ 

WRITE(6»19) { EMH! I»J»1)»I=1»10) » (TRH( I » J » 1 ) * I = 1 , 1 0 ) 

19 EukmaT M0(2XF8,2 ) //) 

12 CONTINUE 

'TTU^s ^ 

END 

FDR I RAT) NI5Kf.K :: 

SUBROUTINE- 0 1 ELCO ( FO » TG » S ♦ El *E2 ) 

PlN2 = Ar71N( 1. )/2b. 

DAM=30./F0 r ■ ’ 

SWN = S/SB.4b ^ ^ 

TC=TG-273. 

F5=' 87 ;"8^ 1 b , 3* 5 WN-0 . 3 A ’i* T C 

DMS = 3i38-0. 1 l*TC+0i0ol47*TC-»TC + 0.0173*TC*SWN -0.52*SWN 

mrG=bT»'SWN-(-0 . 1 2* TC^SWN + O'TTJ'A^'C 

DN=1.+PIN2*(DMS/DAM)**0.98 

UD = 2..*D7J-1 . + ( 0M5/DAM)#iFl.y6 ^ 

El=4.8+( ES-4.8 )#DN/DD 

= l8. *T HKO/l-fJ + l fcS-A»a ) * t DMS/I.)Am ) **0.98/DD 
RETURN 

F-NTi ^ 

S EXECUTE 

T TTTxrrs — io»20ono» »1 odod ^ ^ 

S ENDJOR 






A*2. A computer program for computing the surface emissivity^ due 

to small scale differential scattering coefficient < ( 6 , . 

¥ 


OPTION FORTRAN»FRrNT/15/ 


DIMFNSION A(61 »61 ) .P(61 ) .E(61 ) »REFH(9) »REFV ( 9 ) ♦ TAH ( 9 ) ♦TAV(9) 


l I V ) 

COMMON /ZZ/VCO»THEI .THES»PHIS 




DATA (C0( I) ♦ 1=1 *A ) /4. ,6. ,7. *8./ 


PI 1 0 = PI /TO. 


t(XP = 7*KP-l 




El =49. 


DO 4 11=1*1 

CK T i^CO I in 

El 2=(E1-1 . )*»2+E2*E2 
DO ^ C=TT8 ^ 

THEI=PM 0*FL0AT (K-1 ) 

SI I =SI N I I HR 1 ) ^ — 

CT'l=rOS f THFl ) 

CALL RFFCO r r i -I * 0 1 ..! *pRl*R2I *B! *01 ) 



At(MH = 0.0 


PH1S=F (M) 


Av=n,n 

vro =n . 1 7 

r>o 17 17 = 1*7 





rOMT INIJF- 

IN 

1 ? vro=-vco 

SliMV -SUMV+ .?5«-AV>f S I N ( THES > *A ( L ) 

3 fONTTNU P 

REFV( < ) =.SUMV*.06?5/PI *FCTR1 

Rp-mi — - K - .n6?3/Pi»FC T Ri 

tah ( k: ) =tg*rffh( k > 


TAV( <)=Tri»RFFV( <) 

? CONTINUE 

WRITF(6»5) REFH,RFFV»TAH,TAV 
5 F0RV|ATnH0»9(?)^En.4>//) 

4 (ONT INUF 

STOP 



SU^Rr>lJTTMF REFEOf 

/Yx/EU=T' 

PhI2=ATAN(F?/E1) __ _ _ 

■"'g'= m * '( f-: 1 E 1 + e 2 = 9 ’) * «•", 5 

FU = R^--^- . _ 

' ■ ■pi.TTo=V5»ATANYF ' 

P = ((-. ' J * E + E ? « =■ 2 ) * * . 2 5 

P = p«C0'S ( PrT ri'~) ' ' 

0 s ? <* s I M _< p H I n ) _ _ _ 

PRP2=ljai; + 2'. « J»‘5«r;n5 ( PH T P ) + P2 

■nRL2rP» :.r+2 ; * 9a ’«cn." ( ph i j-pm I n )■; Bis , 

Rj-;.Ty Pjv 

gmd 

^ FO - R - T - f? - AA | \ i DECK 

SUBROUT I NF ANGL F(A*S»T,U»V»W»X,Y»Z) 

fOM f ^Ofl— /■Z " Z/Vn3 - . 1‘ HF I »THFS »PHI S 

B=SORT ( 1 ,+VCO*VCO ) 

S- ( rO S T T HF t ) H VCO* - S - I ' N ' < "T t- t E1 » *rOS" t A > - >7P 

T = S0RT ( 1 ,-S<^S ) 

U= - SIN ( A ) x.SIN ( T H E T") /T 

V=SORT ( 1 ,-U*U) 

W- ( rOG ■^ T M FS) “ VfO^<S!^l < T M FS ■ ) ' »COS ■( A ■ PHIS) )/B 

X = SORT ( 1 ,-W*W ) 

¥ = ^ -f- N -(- P l 4IS"A) « SlN( - T H FS -) -/ X ‘ ^ ^ 

Z=SORT ( 1 ,-Y*Y ) 

RimrR-N ^ 

END . - . 
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SlJf’ROUTrJE CTE'‘F'^(A,P,'=iPl.K,K'^|fCT:?1,<Y) 
T- 1 Vh 'o S n i'i' ' A ( .< Y V Y V'V ^ 

r = K-i 

DO 1 J=?,KX»? 

.7 ‘ A( J,‘<)=-^2.0 

t A(j,:<-1 ) = 1 n ; 

Do 2 J = ^;<x» ^ . 

. ) = 3.'’ t n 

2 A < J, K-1 ) =-f,4 , 0 
_ A(1 ,K)=7,0 - 

A(KX,K)a7.n 

A ( 1 , K-1 ) =- 3 2 . ■) 

A(i^X>K-i) = =-2.n 

Df) n Isl.KX 

"do li J = 3,<.? 

_L = i<.-J . 

A<i,L.fi)sp; j»Af 1 ,K) 

If g. EO.Q) "D TO 11 

A < T » L ) = A ( I , i< - 1 ) 

FCTRt s4 , / ( 900 , QiiC*C) 

_0 = M/C 

r(i) = o.o 


no i> i = 2 iK 
=5 p{ i )=D( 1-1 )*n 
Do ^ I = ^ * K '< 

6 pi. iA.gRji 

Jl = 2 

_ DO 1? IsKl.KY 
; Do 14 Jrl,'<x 

14 A ( J , I > = 4 ( J , J 1 ) 
tS 

p(';Tr-':'o ' ' 

f- r 

ETFTTTTT 

f' LIMITS 10*48000 *» InOOO 

T PNOJOR 

»»» f o F 
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A"3. A computer program for compuHng fhe modified Fresnel reflec^ion 
coefficienf. 


IDI-NT 7 *5 7 , WWU SHTH f S F’ FTO ^ ^ ^ 

$ OPTION FORTRAN, FRONT/l ^>/ 

T M)R I~ R7 T N ~ I 7IT >1CI^ ^ ^ ^ 

DIMFNSION F(5 ) ,SH(9) »5H1 (9) , SHJ { 9 ) » T ( 9 , 9 ) , SV (9 ) , S V I ( 9 ) ,SVJ(9) 

rOMMON " ■'^r r'TT; , CP7 , .9 P7 »'F I TK ^ “T" 

EXTFRNAL FHl , FH? , F VI » F FX P , FV? 

Data {F(J),jsi,5)/7 , < 37 , ^ ^ ~ 

El =54 • 3 . 

E7=T8T5 ^ -^- 7 ^ — ^ 

PI=ATAN( 1. )#4. 

on 5 ! = 1 , ’4 ^ ^ 

THI = (FI.OAT ( n +5.5 ) «Pl / 1 R. 

FT 7 = c, I M ( I M ! ) ^ ^ ^ — — 

SN7=2.*5 IN ( THI ) 

FHM I ) =M.O ^ ^ ^ 

.SHJ(I)=0.r 

5VM 1 )=n.n ^ ^ 

svjm=o.n 

DO 6 T =T ' , '5~ ^ ^ ^ ^ ^ ^ 

PHI=FL0AT( J-1 )*PI/16« • 

7f'=.FN7*( OS I PHI ) — 

FF=FN?*5IN(PHI ) 

CP7=( nS( PHJ )**P ^ ^ ^ — 

FP7 = .FTNf 7.*PHI ) 

ATT- = -nTD : ^ : 

AT2=0.0 

AT 3 = 070 ^ ^ 

AT4=0.0 

fts n.O : 

DO 7 L=1 »4 

B gA+l. - Q- ^ : — — ^ ^ ^ 

CALL TR0MB{8,A,R»FH1*FEXP,T,5»9) 

r i= T (5 - , 5-1 ^ ^ ^ ^ — : 

CALL TROMPf8,A,B,FH2»FEXP»T,5,9) 

r ;;s T (5,5 > ^ 

CALL TROMBt 8,A,R,FV1 »FEXP,T ,5,9) 

T T = I ( 5 V 5-) — 

CALL TR0MR(8,A,R,FV2»FEXP,T,5,9) 

— 14 = I (5,5) ^ ^ 

WRITE(6,2) A,B,T1 ,T2»THI ,PHI ,T3,T4 

' A 'T l = A T l ' -f T n ^ ^ 

AT2=AT2+T2 

AT 3 = A T ^ ^D r— T— ^ r— ^ ^ — — — 

AT4=AT4+T4 

r-- A = R : — 

SHI ( I )=5HI ( I )+ATl*F(J)/90. 

SH j r I ) ' = -S Hjn ) + a t7 *f ( j r / 9Tr; 

SVI { I ) =SVI ( I ) +AT3*F( J) /90. 

svji ir=svj'n )+-ai4 »f(J)'/9u. 

6 CONTINUE 
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S H fT t gFyP f-. ' S T ? ) / PT*r SH T f l T V ) ^ i 

SV(I)sEXP(-ST?)/PI«(SVl ( n+SVJ( I ) ) 

C ON T INUE ^ 

WRITE(6»2) SH»SHI tSHUtSV^SVI *SVJ 

2- F O RMA T 19 t 2 XE li.» I / /) ^ ^ 

STOP 

PTTO ^ — ^ ^ : — 

$ FORTRAN NDECK 

FO NC TTO N FEXP (X ) ^ ^ ^ 

FFXP=PXP (-X*X )*X 

RETURN 
. END 

3 FORTRAN NDE'Ck 

FUNCTION FHKX) 

COMMON St *S.S .CP2'rSP7VEr*E2 .5T2 

. S2=SORT(2.) 

rs-scTx ^ ^ — 

P=Y/S2 

Z= 5 S*X ^ ^ 

. . 0=Z/S2 

r OSHPg ( EXP (P )+ EXP (-P Vr/2. 

EP=EXP(Y)-*-FXP(-Z)+2i*COSHP *EXP(-0) 

T F - ( X if Evl .0 ) P-SQR T ( 1 .-X«X ) ^ 

IF(X.GT.l.O) R=0.0 

E 1 ^-REAL ( C5QR T( F 1-5 T 2 f 0 ,0 > 1 ; 0 V» E2' r i 

CR=REAL(CSORT(E1-X«X'(0.0*1.0)*E2 ) ) 

FM l-:-EP * <B-»-ER^CR ) ^ ■' — — 

RETURN 

Ett© ^ ^ 

S FORTRAN NDECK 

— Ftt NC TTO N F H2 fX 7 ^ 

COMPLEX RtC 

CO MM O N - Sr«SS»CP2*S P 2» F l « E 2*r - . T 2 — — ' 

52=SORT(2.) 

Y=- S C*X 

z=ss*x 

p -^ r rs ?. ^ 

Q=Z/S2 

C 05HP = (EX P(P ) -Fg >t Pf^P - > -r /2. ^ ^ 

SINHP= (EXP(P)-EXP(-P) )/2. 

EP ^ - XP ( Y ) » ( r. -CP 2-) -»-CP 2- »EXP ( --Z H-EXP ( -Q ) * ( C OSH P +S P2 *5 I N Hp-^ 

TF(X.LE.l.n) R=SORT(l .-X*X ) 

IT ( X . G T VT . 0 ) R = rcr . -O » IVO 'I »5 Q R T t X » X-1 « 0 ) ^ : 

r--CSQRT ( E1-X»X-(0.0*1.0)»E2 ) 

EH 2- RF ALTTr^n 'T7 T X- *X+R»C I )*X*X* FP 

. return 

— — FNO ^ : 


A 
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5 FORTRAN NDFCK 

F UN CTI ON FVTf X) ^ ; 

COMPLFX n*C*BCS»CMB 

COMMON BC ♦ S 5~ » C P ? ♦ 5 P ; ? * E 1 > FT7 » ST2 ^ 1 ^ 

SP=.SORT(?.) 

: — — : — ^ ^ — 

ECS=(El-(0.0»1.0)*E2>/( (El+1 . ) *CT2-1 •- ( 0.0 ♦ 1 . 0 ) *E2 *CT2 ) 

r= ;t;r *x ^ ^ 

p=Y/s7 ‘ ; 

Z=SS*X ^ — 

0=Z/S2 

C051HP= lhXP(P)4-hXP(-P) 1/2. — ^ : : 

EP=EXP(Y)+EXP(-Z)+2.*C0SHP *EXP(-0I 

I F rx.LE. 1 .u T- H=baR 1 1 r ' .-x*ri ^ ^ ; 

IF(X.GT.1.0> R=(0.6»1.0)*SORT(X*X-1.0) 

c= cs OR 1 1 E 1 - X » X - ( Tj; 0 , r;' 0 ) *e, 2 ) r 

CMB=ECS*(CSORT (E1-ST2-(0.0»1 .0)*E2)-B*C»ST2*(C-B) / (X*X+B*0 ) 

F V r= R E A L ■( B FCMR -CT# E P ^ ^ 

RETURN 

END ^ 

■5 FORTRAN N0EC< 

FUNt r [ON FV2 t X> 

COMPLFX B»C»ECS.CMB»SEL»SEO.CFV 

~ COMMON SC»SS7cP2»BP2»E1 *E2.ST2 “ 

S2=50RT(2.) 

STl’^OWrtSTJI 

CT2=1.~ST2 

5FL^T2/(E1-(0.0.1.0^*E2) 

ECS=CT2-BEL 



z=s.s*x 

p=V/S^ 

0=Z/S2 

— cosH P = (FX Prrn -t- E X P T- P T i /2 ~ ; 

SINHP= (EXP(P)-EXP(-P) )/2. 

C P l=S0 R['t CP2)* (F y P ( Y)+B2» ETP (- Q >*s r NK P ) 

SP1=S0RT ( I .-CP?)*(EXP(-Z )+S2*EXP(-0)»C0SHP) 

EP =E XP ( Y ) *C P 2 + I 1 .-CP2 1 * E X P~r= Z ) + E XP [ -0 ) * ( C0bHP-b P 2*b 1 NHP > 

IFfX.LE.1.0) B=S0RT(1.-X*X) 

1 F ;g t . r. o r- B ' =f o . o » i ; o ) *sqr t ( x*x~i .0 1 ^ 

C=CFQRT (E1-X-*X-(0.0»1.0)*E2) 

SbU=( SUR T ( E r^>T2- 1 0.n"i 1 . 0 ) *t2 ) 

CMR = X/(X*X^-B*C) 

C F V= ' CMB*( 1 .-SEU*X*E P -2.*S [r *&E0* ' lLPi~SPl ) ) /bCS-- 

FV?=RFAL (CFV) 

RCTTTITN 

END 
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F O RTR A T r -NrrFCK : 

SUBROUTINE TROMR ( NMAX , A » R » F » G,T»JMAX»NRC) - 

D T M E NS T UN T( N R C»N R C 'I : 

H=B-A 

T (1 « 1 ) = 1 F ( A '1 »G f A ) -t- F ( R I «GT R-) ' 

DO 2 N=1»NMAX . 

n N+1 ♦ 1 ) =0.0 77~ ~~ ~ r ~ ^ 

FR=H/2 ,n*»N 

rwAT^T^nFTrn ' 

DO 1 T = 1 ,TMAX.2 - 

1 T ( N-H ,rr = T (N-H .1 )+F(FLOA T t ! ) * FR +A )»G( F LOA T ( I ) * F R-<-A ) 

2 T ( N + 1 , 1 ) =T (N. 1 ) /2.0+H*T ( N+1 » ] ) /2.0**N . . 

DO 3 J = 2.JMAX T : — ^ ^ ^ ^ — : 

NXMjP2=NMA5<-J + 2 

FOR-JWT-4;0»*r J-1 I — — — 

DO 3 N = 1 .NXMJP2 

3 TI N.J)^ (F 0RJM1* T TN+1-»J-1 ) - T (N.J-1 ) >/(F0RJMl-1.0 ) 

RETURN . .... 

ETTD ^ ^ 

f> . ■ FXFCUTF .... 

■'* L T M T r .s ro'nTTnTnT'iT i a n o n ^ ■ 

$ ENDJOR . . - • < 

--» ^ - E t) F ^ ^ ^ — : 




A”4. A computer program for computing the two~scale rough surface 
emissivity 

If DFNT 75TT7FTr’TTFVr~i7T) ; ^ 

$. . OPTION F0RTRAN»FRCMT/T5/ 

T FORTRAN NHFCK ^ ^ 

DIMFNSION A( 61» 61 ) ♦P(061 ) *E (061 ) »TATM(U9) »ATZ(09) »TREF( 10»3»9') 

1 *XFS1 1 061 ) S:^ CD'S! nTFS’ll 06 I ) »(.E"SAT06 1 ) »CF6b (061 ) ♦CFS6r06ri 

2*SFN1 (.061 ) » SEN? (061 ) *SEN3(U61 ) ♦SENA (061 ) ♦SENS (061 ) ♦SEN6(061) 

.5*VEL C3 ) ♦FMISTiUr3»b ) » I AK 5TTU', 3 ♦ b ) ♦ I MES ( iO ♦ 3 ♦ b ) ♦ TSKV'I 09 t 

AfCES7( 61 ) ,CES8 ( 061 ) ♦CES9(061 ) ♦CES0( 061 ) ♦SENV ( 061 ) ♦SENS ( 061 ) 

5*Sir.N0( 61 ) ♦seNv ( 0611 ♦AP I 61 ♦S] ) ♦PPi 61 ) ♦eP ( 6T) ^ 

COMMON PUPI360»F1 ♦EZ^CES^SEN^EUI ♦E12^ES^EC^PE^0E^TM^S20 • ' 

i>Ai A ( VFL(K.) ♦< = !♦? )/l6.0"20.00/ r 

PT = ATANd. )*A. 

PT 1() = '->1 / IS.O ^ ^ i 

PI360-?,0*PI 

PI 0() = p J /^,() ^ ^ 

PIA6=PI90/2. 

PT2T=mr577-; ^ ^ ^ 

SL0P1 = ,?5 

SrOP2 = . 12b ^ ^ 

SL0P3=.0625 . ' 

sropo= . 1 H / b ^ 

TG=291. 

tM=1 ,1?*TG-60. ^ : 

e;?n=, 01921 

— fl=bH..6 ^ ^ 

E2=36.6 

hin= fFl-1 . )t^W2 + F'2^E2 ^ 

FI 2=F1 *F1 +F?*F2 

PHIO=,b*ATANrF;'2 /rT ) ^ 

ES=SORT(F12) 

Pir= SOR TIES! ^ 

PF=RO*CO.S( PH 10 I ■ 

^ ^ 

OF=RO».STM(PHIO) 

irr=Fn — ^ ^ 

,<YP = 61 

KTTTl 

<P = 20 

K5(=2*K-] ^ ^ ^ ^ 

<XP=2»<P-1 

CALL C0F F FS(A ♦? ♦ E * P dK ♦KX ♦KC. I RNK.Y ) ^ 

CALL COFFFS(AP,PP^EP»PI ♦KP',KXP^FCTR2t(CYP) 

CALL SORT IP^OlVPiyO^SLUPdPl ♦CtSl ♦bENl ♦KX ♦KY ) ^ 

CALL SORT(P^PI A6^PI9n»SLOP? ,9! ♦CES2 ♦SEN? ♦KX ♦KY ) 
e ALl mjk I j PP ♦ 0.0 ♦P I A2 ♦bLlJp^; ♦K 1 ♦(.tbi ♦i.t: N 31'K.XP'^ KYP ) 

CALL S0RT(P^PI90^PI ♦SL0P1»PI ♦ CESS ♦SENS ♦ KX ♦ KY ) 

< A L L .SORT ! T ♦"“n'. 0 ♦ P I 2 2 ♦ SL OP '3 ♦ P I ♦ C b S9 ♦' S h N 9 ; K X ♦"K Y 1 

CALL , SORT(F^PI22^PI9n^SLOPO^PI ♦ CE SO ♦ SENO ♦ KX ♦< Y ) 

nrr~9rrnrsT79 ^ 

THFTA1= FL0AT(K-1 )«pll0 
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r ^-T O'; CTMt- T'A'I ) 

<,FN = r,IN( THFTAl ) 

XT=THFTA1-Pr7^ 

IF (Xl.LT.n.O) Xl=0.n 

— x7 = thf lAi+pi;?? ■: 

IF ( X? .GT.PIPn ) X? = Pl9n 

SLI)PA = X 1 /PI 'Hbn ' ^ 

SLOP')=(P I9D-X? ) /Pn6n 
57(-OP6= ( X7-XI ) / P I '^60 

rALI, 90RT(F ,o.n»XI ♦SLOPA»PI ♦CFSA»SFN4»KX »KY ) 

CALL SORKF- »X?»PI90»SL0P5»PI .CFS5»SEN5»KX ,KY ) 

CALL SORT(E»Xl»X2»SLOP6»PI »CES6 »SEN6 »KX »KY ) 

CA'LL" £0RT(EFrXT,x'2, — 5L0P6';P I ,CEs7, SeN 7 ,KxP ,KYP ) 

00 90 N=I »? 

DO 90 ^ ^ ^ I 

C0 = . 004+0. 78E-^*\/El(KK ) 

FF2 = .ALOPl*.<!irOP4 ^ 

Pl=PINT(Pil*SLOP4»CES4.SEN4,CESl»SENl,CO»N, FF2.A »KX »KY ) 
FF2=SL0P1 *SL0P8 

■ P2=PINT(P22»SLOP5.CES5,SEN5,CESl.SENl»CO»N. FF2»A ,KX ,KY ) 
FT?=FL0P2»SL0P6 

' ■P8=PINT(P33»SL0P6»CES6»SE'N6.CES2»SEN2»C0,N» FF2»A,KX»KY) 

FF7'=SL0P2^(-!^D:7P1?) 

P4=PINT ( P44»5LOP6.CES7,SEN7,CES3*SEN3»CO»N, FF2 » AP , KXP » K YP ) 

FF2 = ?,LOPl^-.3i:OP3 

P5=PINT ( P55.SL0P3»CES9»5EN9,CES8»SEN8,C0,N.FF2»A»KX»KY ) 
FF2=SLOP1*SLOPO ' 

P6=PINT (P66»SL0P0 »CES0.5EN0.CES8»SEN8,C0,N»FF2»A,K^*KY ) 
WRlTE(6»7) PI ,P2.P3*P4,P5»P6 
7 E0RWAT(6(5XEn.4)//) 

T1NT= { P] +P2+P34P5+P6' /PI360*FCTR1+P4*FCTR2/PI 360 
TAREK ,N ,<K ) = n .-T INI ) *TG 

TKhF ( K, ,N *K.<. ) = I ( F 1 1 +P2 2 +P33 +P66 + P66 ) .*FC T R 1 +P44*FC TR 2 ) /PI 360 
1 / ( TM* n ,-FXP (-S20/CES ) ) ) 

IMESI <,M ,)0<) = AI 2( <)*n ARbCiCTN'.KK ) + [ RLE ( K.,'N,KK J ) + TATM{K ) : 

90 EMIS(<.M»KK)=T.-TINT 

00" I 5 N= ! .2 ^ 

DO 14 L = 1 ♦ 2 

1 4 w irTTr rs . 5 ) eei ^ S' i k . n ) . k= t' V9 i. > u ar s ( k » n ,ot< ' = r« 9 ) — ^ 

I * ( TREE (K.N.L).K = I,9)» ( TME5 (K »N»L ) »K = 1 ,9 ) 

S h O 'R' M ' A n I HO ^ ^ TT P ' X ' F T '7 4 ) T71 

FMIEno.N.L )=0.0 

1 PI- H I n.M.L r= 1 .0 

18 rOMTINUF ‘ 

WR T T F I 43T70 r‘- I ( (F7TT3‘tXT7^inTK = til 0).N=1»2)*L = 1»2) ^ 

1 ♦ ( ( ( TREFK.N.L ) »K=1 »10) »N=1 .2 ) »L = 1 »2 ) 

70“ FORWATM 0 1 F BT^ni 

ETOP 

rMn : — : 
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S FORTRAN NDFCK 

STTBROOT TNF~OORT ( F ,Z T»Z? ♦ SLOP »P I tCFS i SETN .ICXVKY ) 

DIMFNFION F(KY I »CFS(KY) *SFN(KY) 

IF (SLOP. FO, 0.0) GO TO 70 

DO I 7 I -I ,KX ■ 

x^rrm-) 

ANG=Z1 + f X+P I ) »SL0P ' 

rF"‘;( I ) =ro«rrANG) — ^ 

17 SFN( I ) =G IN { ANG ) 

VO "RF TURN 

END 

" j.VRnuT i v = .., a , p , r- , = I , i< , K < . >='n tri , < y ) 

n r H E V s 1 D 'I ^\ ( < Y . < Y ) , ^ ( K Y ) , r ( i< Y ) 

M = 3 I 

n = < - r 

T 3 1 J = 2 . < X , ? 

A ( J • O = - 1 2 . 3 
1 M J» <-U =1 30n . n 
3 3_2 J = 3,<X',2 
A ( J . ■<) = 1 't . 3 
A ( j . < - 1 ) = -. 5 1 . n 
A (1 . Y ) = 7 . [/ 

A ( <X,. Y ) = 7 . :i 
A ( 1 » > - 1 ) = - 3 2 . 0 
A{<x,K-n = -3^,;i 
33 11 I=li<X 
^3 11 J = 3 » < . 2 
L = Y-J 

A ( I , L 1 ) = 2 . 3 "• A ( I , K ) 

I" (L. =3. 0 ). 30 T3 11 
A ( f , L ) = A ( I , < - 1 ) 

ll" C3'lTt vj j= " ■ . 

r 3 T R 1 : 4 . T r- H a -I / ,( , 9 n P . n > C C ) 

' 3 = -l/C 

>= (1 ) = 0 . 0 - ' , 

'^ ( U = -•) 

P3 5 1=2, < 

- ( ! ) =E { 1 -1 ) +3 , . . 

? I )=P( 1 -1) ^-3 
X'] = k + 1 ‘ 

P36I=<1,<X 
A -•> ( T ) = P ( I - 1 ) <■ 3 
,)1 = 2 

33 I'-i I ='<1 , <X 
3 3 I 4 J : 1. , < >< 

'r'4 . 1 , 1 )“= A ( j , j 1 ) 

F (!)=-=( .11 ) ’ 

lT,llsJl + l, 

PE TUP N 
c;nj'i 
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■n FD R T RAN ND F CK 

FUNCTION PINT ( Q , SLOP » CT H? ♦ STH2 » CPH? « SPH2 * CO » N » F2 * A »KX » KY ) 

DTMFNSION ATTYtKYI »CTH2(KY) ,S’TH2(K;Y) ,CPTI2{KY J »SPH2 (KY) 

COMMON PI *PI360*E1 »E2 *CES»SEN.EUI » E 12 » ES » EC ♦ P E »OE ♦ TNI » S20 

Q - 

PINT=0.0 

IT — (SLOP.FO.OVn ) GO TU 20 

IF(CO.LT..016) TTT=1*754*BS(63.0) 

IFICU.GE. .OT&T I T ] = 1 • 754*B5 ( 95 . CD 

DR !=ES+2.*PE*CrS+CES*CES 
DL"T=ES+2.*E(.*CE5+CES*CES*E12 
- DO 17 I =1 ,KX 

Cl I=CIH2( 1 I 

TSKY = TM*n .-EXPI-520/CT1 ) ) 

~ SI l=STH2 I I ) 

V3=CE5+rTl 

VP=2,»V3**2*C0 ^ 

VZ2=V3*V3 

DRS=hS + 2.*P'E^I 1 + L M^C I J 

DLS=ES+2.*EC*CT1+CT1*CT1*E12 

DO” 17 M=1 ,KX 

SP1=SPH? (M ) 

CPT=CPH2 ( M ) 

V1=5T1*SP1 

*rPl-StN 

VXY = V1*\/I +V2*V2 

V^»=M I •*(.'ES*CP I +SEN-»(.Ti 

Rl =V4*V4+V1*V1 

U2 = | .0+CES*"Crr-SEN*5! l*rP] 

UU*U2*U2 

U=SIJKI IU2/2.U ) 

IF(N.F0.3) GO TO 14 

klT2 = l,-4.*PE*U/(E5 + 2»*Pt*y+U'-«'lJ5 

RV2 = 1 .-4,»EC«U/( ES + 2«»EC*U+U*U*E1 2) 

RH:J = KH2^^( r.-874U*TTT-*l-KCH(U)") 

RV? = RV2^M .-8.»U*TTT*FRCV(U) )• 

— 14 (TmiNur 

IF(N.>^0,1 ) FH2=(RH2*V4*V4 + RV2*V1*V1 )/Rl 
I F'IN.FO.? ) FH? = (Rv?^V4*V4 + RH2*V1*V1T7T^1 
IFIN.FO.3) FH2=l.n 

GAMA=4.«UU* FH2*EXP (-VXY/VP ) /(CES*VP*VZ2 ) 

SPS2-SPl->t,SPl 

CPS2=CP1*CP1 

OT = SEN**2-2. * SEN*ST 1 *CP 1 +ST1*»2 

p«K=4.^TTT^EXP(-or )4cFs*rn**2*Eui#4‘; 

TF(N.FQ,2) go TO IS 
I I N*PfJ. 4 ) GO TO 1 A 

GAMa = GAMA+PMR« (rpt:,?/nRS + .SPS2*F5/DL.S ) /DR I 

fTO" 10 I « ^ 

IS VNl = ( n *SFN*ST 1 -CPI * ( PF*PF-OE*Q) ) )**2 

VN/'= ( 2 .*CP1 •»PE«T5F-F2*SFM*1^'TT )^<■*2 ! 

GAMa = GAMA+PMR» ( (VNl -*VN? ) /Dl. S + SPS? *F S/DRS ) / DL I 
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T 


lA CONTINtlF 


i « l UIN 1 1 

PT INT = CiAMA*STl I ) 


pint=ptnt+ptint 

■ 



17 continue 

?o rfturn 


; V 


FN13 

* FORTRAN NOFCK 


. - - 

• — 

FliNOTTON R<^(X) 

RS=A.n‘^F-^*n ,?A7+0.r.?68*X+6.03E-5*X*X )**? 

nr- mm — 

<■ 

. ' ■ ■ 



r\ 1 1 1 rv 

END * . 

t FORTRAN 
FUNCTION 

NOF.CK 
FRCV( X > 

. , ' J 

IF (X*(il • J 
IFTX.LE.l 

• 0 •0r,X,LT.1,E'' 
.O.AND.X.GT.0.5 : 

-b) FKCV = 0.U “■ ~ 

1 FRCV=1./(8.57*X-4.4*X*X-2.243) 7 

IF CX.LE.O 
I F ( X . L E . 0 

. 5 . AnD.X. ij 1 . u . 1 , 
.173) FRCV=l./(] 

Tn FRCV=1. / ( 0 . 8 +X+ 2 .Bi*x*X-U. 14 ) 

LI .1*X-55.3*X*X-0. 184 ) 

RE TURN 
END 


t" OR 1 K A N 
FUNCTION 

NOtCX ^ ^ ; 

FRCH(X) 

IF(X.M.l.O.UR.X.LT.l.E-) 
FRCH=0. 7?4A-)f x->^*7 ,5+0.294/ 

i) hO JO i ^ 

00 TO 2 
1 FRCH=n.n 



FXFCtiTF 
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A-5. A computer program for computing two scale backscattering cross" 
section, 10 ) . 


*. . . IDFNT WU SHIH TSENG 

t OPTION FORTRAN ♦FRCNT/15/ 

FORTRAN NDf-CK 

OTMFNSION CT(3).RMS(S),ARH(9.3)*ARV(9,3) »STGR (9,3)»TCAH(9,3), 

I K, AV I V » 3 )■ , A IH ( 9 , 3 r, ATV ( 9 , 3 ) 

~DAtA ( RMS ( J ) »CT ( J ) »J=1 »3 ) /,10 , . 1 3. . 12 ». 17« . 15».20/ 

FT^TTBnr? ^ ' 

jr?=34.9 

CALL Rb'K’UI 1. ^ 

■"RFs=r.-4.*ZA/Z2 

2 - 1 » 9 ^ 

THFI=PI10»FL0AT (K-1 ) 

V » M- r F?yt2. - -C5 2 ) 

VN2= (F1*(2.-CS2)-1 .+rS2 )**2 

DO ? M-1 *3 ^ : 

CMF=16.*CT(M)*CSA*FXP(-A.*(1.-CS2 ) )*CT (M) 

G 2'= RM? ( M t ■ . 

Vr0=2,*G2*CS2 • 

STGR rr, ,M-) r:RFS/ ( venues?, l *F.XP l - ( 1 . -CS2 )7 VCO I * (1 , -4T»T:T ( M ) **2 ) 

■'2 CONTINUF 

P O A J = T > 3 

VN=RM'5( Jt/S0RT(2. ) 

Vn-=5iOR T ( 1 .+RMS( J) *^»2 > ^ 

DO A < = 1 » 9 

T tt g f ^I^ MF- tOAT (K-1 ) 

RH=0.0 

R VgC Ti T) 

FA=1.0 

DO 1 2 T 2= lt ? 

CXP=(C05(THFI )+FA*SlN(THFI )*VN)/VD 

C< ;4 = rF ? *CS2 

rOMF=«,*CT ( J> «-*2*rSA*FXP(-4.»n ,-CS2 ) ) 

R H = R TTT-rnM FyR H2 

RV=RV+COMF*{ VN1 + VN2 >*E1 ?/DRL2**2 

12 FA = -1. 0 

ARK(K»J)=RH 

ARV t lCO V=RV ^ ^ 

ATH{K* J)=ARH{K» J)+STGR(K.J) 

A T V (R » j -c =ARV ( K » j r rs Trr R ric ^ J T 

4 CONTINUE 

WT T F ( 6 - » -? 4 )- S T GR FA T H »- a T V > ARH t AR V ^ T C AH « T CA V 

74 FORMAT(9{?Xrn ,4)//) 

^ F X F CL ITP ^ 

^ LIMITS io.?nnoo» »lnnon 

■T F N nO OR 

**-»FOF 


131 



A“6. A compufer program for compufing the clear sky and cloudy sky 
apparent temperatures. 


■? TftTKlT 75?3»W^»WU SHIH TSFNG 

$ OPTION FORTRAN,FRCNT/15/ 

« FORTRAN NOFCK 

S INFOOF IBMF 

niMFNSION FMtsno ) »RFrOMO,2»4 ) » TMFS (1 0 . 2 » 4 ) »C ( 100 ) » 

irPT n 0 ) ,rTP( J 0 ) «rLnM( 1 0 ) » AZ ( 300 ) ♦ TZ ( 300 ) ,z 1 300 ) » ATZ ( 1 0 ) *TSKY(1 0) ♦ 
~ ? TATM(loT,Xrmno ) .CSPF( IO) ■ 

PI =ATAN(1 . )Hc-4. 

pno=pi/i8. ^ ^ ^ 

FR=13.9 



H=7.0 

P M = rZ 30Z6 ^ 

AIC=-4. 374*1. E-4 

RTC=-I . 7-67*'l,T-'2 

RFAr>(^,3) ( ( ( FMIS(K»N,KIC ) ,K = 1 , 10 ) »N = 1 ,2 ) ,KK = 1 »4) 

r» rfTRTrOnCVN,XKtTK = l , 10) ,N = 1,2) ,KK=l i4) ^ ^ 

3 FORMAT (T0(F8.4) ) 

CA^ - t r ATM F FTT G , r ;5? T Sn. T ir. 9 ,H > TA TM , A T~ Z ,S 20 , T SKY , AZ , T Z Z ) 

DO 90 K= 1,10 

Dt3 90~ N - 1 F2 

DO 90 KK=1,4 

rWf'St IC^ N ', KF ) = ATZ nc) *■( TG* EM 1 s ( 1C , N , KK ) +RFCO ( 1C , N , IOC) ♦ I SiCT DO^TT AT M FK ) 

90 CONTINUE 

WR - ! T E ( 6,59r 1 S , R F CO , T M E 5 

55 F0RMAT(10(2XE11.4)//> 

[JO 8 - FT = t~ , ' T0 ^ 

READI5,123) CBT(M) ,CTP(M) ,CLDM(M1 ,CSPF(M) 

— 173“format( 4 1 f ro .t, ) 1 

WRITE(6,18) CRT(M) ,CTP(M) ,CLDM{M) ,CSPF(M) 

rR-FORMlV 1 ■[ I RO ,'4' I TOXl KTU"? 477 H ^ 

ICB = (CBT {M)+50. ) /lOO.+l .0 

K 1 = (7 T p I NT) --50“} 7Tmr;"n ro ^ ^ — 

CM=CLDM(M) 

IF (CM.LT70T0DT1 GO '10 16 ' 

DO 11 I=<B,KT 

Cn )=0'.0 ^ 

U=TZ ( 1 )-273. 

r ETiTj+8 . r. n". if; (T'n“E.cr ro' T7 

cif =PN*FR**1 ,95*EXP (-6.866*( 1 . + .0045*U) ) 

ormr-77> ^ 

71 CK=PN*FR**1 .006*EXP(-8.261*( l.+BIC*U+AIC*U*U) ) 

26"rn r=cM*r<*t,*BPFiM) 

WRITF(6,2 ) S( I ) ,T2 1 I ) ,C( I ) »AZ ( I ) 

— ^ 2 TORM'AI 1 I OX, 4 1 lUX,frT2?4) ) ^ ^ 

Azn )=A? ( I )+c;i I ) . 

n -' COTI Tl NUE 

16 CONTINUE 

C I H I S I S T O ~ CAL ~ CUL ATE~rR AN5M lTT ANC E ,, T A T M, I Sia ^ 

JJ = 10. *H 
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no 70 j=i » JJ 
AT'- AT+(Vri»AZIJ) 

70 CONTINUF 

WRI IF{6,'6fn ) ^ ^ 

603 FORMAT ( 1 HO , 3X ♦ 7HDFGRFES » 0 8X »1 3HT RANSM I TT ANCE » 1 UX » 1 ] HT / ATMOSPHER ♦ 

TT^XT5nT7TXY7T 

DO 30 K=1 »1 0 

THETA! =FLUAT(K-1 ) *PI 10 

- IF(K.EO.IO) THETA1=8.95*PI 10 

CFS = COS ( THETAl ) ^ 

' ATZ(K)=PXP(-1 ./CFS^AT) 

ATTN(K )=10.*AL0G10(ATZ(K) ) 

ST = 0,0 

no hO I =1 ,JJ ^ 

S2 = 0.0 

DO "50 j=r;jj 

S2 = S7+0,1*AZ(J) 

“Tn OUN ( INUb ^ 

S3=FXP (-1 ,/CES*S7 ) 

.S1 = ST +0.1*S3*TZ( I )*AZn') ^ 

40 CONTTNUF 

T ATM! K ) =TT7-CE.5*'S1 
••■•S4=0,0 

DO /O 1=1,300 ^ 

S5=0,0 

no 80 ij=i,i ^ ~ 

S5 = .S5 + 0,1*AZ( IJ) 

Fro C'OMriNUF 

66=FXP(-1 ,/CE5*S5 ) 

S4=.s4+n, i*sb*izi n*AZ{n 

70 CONTINUF , . 

^ rSKV ( < ) =1 ,/CE.S*b4 
K3=(<-1 >*10 

WRI TF (6.4 ) K3»A TZ' ( I O t T A T MIK ) * T S I CY( i a >A T . T N(K) 

4 FORMAT(5X»I5»4(5X,E16.8)/) 

“TO C O N TI NU F 

no 15 I=»^R,<T 

AZ ( n = AZ ( 1 ),-( ( 1 ) 

15 CONTINUF 

~l THIS' IS TO'CALCUCATE KMTS ♦ T REE » I MES ^ 

WRITF{6,605I 

“Strs F ORM A I ( I HO ♦ 1 b X W H BE GK E, E b V 1 2 X > 1 IH E M M 1 " bb 1 V 1 I V W3 X7 H H I' 7R E 'F L E O E T S 

1 ♦15X»10HT/MEASURED /> 

DO' 44 <=i,vn — : 

DO 44 N=1 »? 

BfT-4'4 KK = 1V4 ^ ^ 

TMFS(K ,N ♦<< ) =ATZ ( K ) * ( TG*FMIS ( K , N , KK ) +RFCO( K'iN , KK ) *TSKY ( K ) )+TATM(K > 

44 HINT INI IF ^ 

WRITF(6,56) TMFS 
56 format ( 1 0C7XF6.7 )"//) 

8 CONTINUF 


END 
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^ FUR l f^&^r ' N I >K K 

$ INCODE IRMF 

■ SUBROUTINE ATMFFf | C, » l)NG"» PG * FR , H , I A I M , A rZTSZD’i I 5K Y » AZ » T Z ♦ Z ) 

OIMFNSION aZ( 300) ♦TZ('iOO) »Z(300) »T ATM ( 1 0 ) » ATZ ( 10 ) ♦ TSK Y ( 10 ) . 

rrsuPTiTn 

RFAL LWO»LWW,I.WO.S,LWWF, 

OAT A -P S F AT ♦ 1 bT: A , n N b h A U A Z Z , t, CB 1) , t. < 1 i b'V(. K K E S ♦ E 'Z 'Z , F 6 0 , m 8 T^FO 77 
1 290,,1. .32.A,.36Z.6,5.23,2.55E-3,22.23t?,60..1 18.75/ 

R T TO = A T ANn ^ ) /A, 3 

JJ=in.*H 

Tr= l . 1 2« T G-5D. 

FRS=FR*FR 

RM22 ' a I APS[ F R-F-22T) ** • ? 

RP22 = (FR+F22)**2 

TT .i -f ,o ~ (A R- 5 (FR - F 6nh* T f '2 

SP60 = (FR+F60)**2 


SM11R = (ARS(FR-F118) )**2 
SPUR = (FR+FllR) **2 
S20 = n.n 
WRITE (A ♦! 2) 

M 2 -format rTHlT28HATTEN. CUth . H NEPER POK~KM "T/7! 

00 10 1=1*300 

XT- s— I 

Zm =(100.*XI-50.)*l.E-3 

rr-rTTrTrTTTTTr-TTTTTTr-zi 

• TZn > = T6-71.*Z( n/11. 

Go -T O " ZA : ^ 

23 TZm =217. 

rON TT Nl TF 

ONZ=ONG*rxP(-2( n /2*205) 

p z -=■ ■ p r ryF~ x R T -a. i u 3 *z t n i ^ 

PA=1013.25*PZ*1 ./76o. 

LW O 1 .5» T ( I S E A/ T Z ( T n**0.5 ) *PZ/PS E A 

LWW=2.58E-3*( l.+.0147*DNZ*TZ { I ) /PA) *PA*( ( 318. /TZ( I ) )**.625 ) 

G0-= 0 « 2 1 »PZ wlW 0* F R5/ T T Z ( 1 > * T Z( I ) *T Z ( r> ) 

GW=1 .*EXP( -644. /TZ( I ) ) *FRS*PA*DNZ* ( 1 . +.0147*DNZ*TZ ( I ) /PA ) / ( TZ ( I ) 

) 125 )' 

LWWS = LWW*LWW 

LW05 ^ ' LWO*LWO • 

A22 = C<22*GW*(1 ./(SM22+LWWS)+1./(SP22+LWWS) ) 

A 60 - = C > C60*G0« r t ~ ; / t- SM60-*^LW05. ) +1 . 0 / ( 5P60 + LW 0 5 ) . O / ' f F R5 ^ +^W0Sr ) 

A118 = c:<1 lR»GO*n .rt/ ( 6M1 18+LWOS ) + 1.0/(SP118+LWOS ) ) 

•A RF S = r< RP‘ S^DNZ* F RSgl-WW7 ( I Z ( 1 )**1 .5 T 

A2(I)= A2?+A60 + An 8 + ARFS 

F 2 n = R2n-TA7. ( n It. T 

IFd.FO.JJ) AT = R20 

W RIT E (6* 11) — I Z( 1 ) *AZ( 1 ) .Z I I ) *620 

10 CONTINUE 

W R I TF ( 6 ;TT ) ~~ H * a t * S/O * I M 

11 format MH0*4( 10XE11*4)/) 

WR IT E ( 6 * 1 3) 

13 format ( T HO *92HTRANG. TATM. TSKY VERSUS ANGLES *) 
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no K = r^Tn 

THF = Fi.nM f K-i ) »PT 1 n, 

r F ( . f- 0 n m — [HF = H.9^*pr rn ^ 

Ft^r = 1 . /ros ( THF ) 

rsK-Yrn-sTwim ;--rxPT-.*5rr*F.?o ) i 

- •• AT7 (K ) -FXP(~AT*SFC ) 

irrsTTm ' ^ 

DO AO I =1 « JJ 

b/=o.n “ 

DO 50 J=I,JJ 

5 0 - 57-52+«T*A Z ( J ) 

AO SI =S1+.1*TZ( I )*AZ( T )*FXP (-5Z*SFC ) 

r ATM rKT=sFr>s i 

SS=0.0 

DO' "6n"'T"JJ» HOn 

SA=n,o 

DO Vo J=JJ»t 
70 SA=SA+. 1*AZ ( J ) 

60 S^ = S54- - .r* T'Z l 1 ) *AZ ( I ) * F X P (-SA*bF(- ) 

TSUP(< ) =SFC*S3 

g:^”TT : - i TyT n ~ 

30 WRITE(6»1A) K3»ATZ(K>»TATM(IC)»TSKY(K) ♦TSUP(K) 
1 A F ORMA T ( 5 X t T 5 » A fl OX F riVAT T Z ' ) 

return 

end ^ 

7 EXECUTE 

5 I.TMITS 10»?0000, ,]'>000 
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A-7. 


A computer program for computing the apparent temperatures of light 
to moderate rain. ... 


I DP NT WU ^ ' 

$ OPTION FORTRAN ,FRCNT/1 5/ 

f. FORTRAN NDFCK 

S INCOOF IPMF - - • 

KfiAL lwo,lww»lwo«;»lww?^'lwv»lw\/s*lam ^ 

OIMFNSrON C(300)»Hno) »RFrO( *FMIS( 10»?»4) »TMFS( 10,2 ) ♦ 

IKriUO) ,AI (.vriui ,F0nn) r,AZ(^00) ,F( 7) ,TZI300 ITatZTiOTI 

2TATMnO) ,TFFYnO) 

romFnrr- ^ ^ ^ 

NA.MFLTST /LIST!/ TG»DNG,PG . ' 

DATA PStA", 1 SEA , DNSt ATCK22 ,CK6aiCKrr8 ,CKRtS , f,22 , F60 , 1 IB / 76u . , 
1 291.,1. ,32.A,362.6,5.23,2.53F-3,22.233,60.,118.7b/ 

GK"183'^8,6t-3 

F183=183.3 

PN'=.23D2fc 

AIC=-4. 374*1. E-4 

RrT=-i. 76 7*r.t-2 ^ ^ ^ 

A890=.0n8- 

ATTT^TnT73 ^ 

AT 39=,0?6 

H89n=T .'^2 ^ 

. Rin= 1.26 

Bl"39=r. I 8 ^ ^ 

PI = ATAN(1.)*4. 

pno =P! /i 8 .o ^ ^ 

El=58.6 

E2^36.6 ^ ^ ^ 

READ(5»13> ( ( (EMIS( I »J,.K1 ,I>1,10) ,J = 1,2) »K = 1»4) 

r» t nRFT:o( r»j»!0 ,i = ivio) ,j=i,2> ♦k.=i»4J 

13 FORMAT ( 10(F8.4) )■ 

REALM 5,L i 3T 1 ) 

WRITP(6,L ISTl ) 

1 hORMA? ( Vh !6;'d 1 

R‘^An(8,10?» fH(J).J=l,3) ■ 

102 F0RMArf3FTn.4') 

RFAD(5»59) NF,NM 

39 FORMAT (2 151 

read (5,239) (F0( J) »J=1 »NF ) 

239 FORMA) ( lOt- K. 2 ) ^ 

DO 8 L=1 ♦! 

WRTT E r(T, 4T"L rFO'dl ^ 

4 format ( 1H0»3X,13HFRF0UENCY FO ( ♦ ( 1 ,2H ) = , F8 . 4 t 3HGHZ / / ) ; 

FR=FOm : — 

C THIS IS TO CALCULATE THE SKY . ABSORPT I ON COEFF. ■; , 

FITS = R*F R ^ I 

S20 = o.n. 

DO 10' T=r»3DT) 

•XT - T • • 

zin =nDn.*xT-5o. )»Kf-3 ' ; 

IF (I. GT. 1 50) GO TO 2 3 
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“ I zi r y= I shA-7i,»/. 1 1 )/i'i . 

IF (TZ( r ).LT.217. ) TZ(I )=217, 

i n 2A ^ r — 

73 TZ( I ) =717. 

zn — roNTTMUT ^ ^ 

r)MZ = nNr,*FXP ( -7 ( T ) /7»7 ) 

PZ = "Pb*hXP( -0. 1 V3*ZM 1 ) 

PA=1013,23*PZ*1 ./76n. 

LV ' O = 1 .3*f (T S F A ' / T Z '<- n ) »*0,5 ) * P Z/ P 5hA ^ ^ 

LWW=2. 58E-3»( 1 . + .0147*DNZ*TZ { I ) /PA) *PA»( ( 318. /TZ ( I) )**.625 ) 

- LWV=3.0*( PZ/ 760. >*(2^0. /TZ( f ) )**0. 5 

GV=nNZ*FRS*LWV»( (290»/TZ( i ) )»*2)/3.0 . 

. . GO = 0.21*PZ#LW0*FRS/ (TZ( I )*TZ( I )*TZ< I ) ) 

GW=1. *EXP( -644. / TZ( I ) )*FRS*P'A*DNZ* ( l. + .0147*DNZ*TZ ( I)/PA)/(TZ(I) 
! ^^*3.125 ) ^ ^ ^ ^ ' 

Lv;ws = LWW*LWW 

LWUS = LW0*LW0, , ^ ^ ^ ^ 

L.K'VF = LWV*LWV 

SM22 = (ARS(FR-F77 ) 

GP22 = (FR+F22)**2 

SM 60 = ( AHS(FR-F60M**2 : 

SP60 = (FR+F60)»*2 

SM118 = (AHSIFR-FllS) )#*2 — : — ^ T: 

' SP118 = (FR+F118 > **2 

bM153 = ( AHb ( FR-h 183 ) ) #*2 ^ ^ T 

SP1.83 = (pR+F^B3) **2 

;■ A22 = t.K22*GW*( 1 .7'rSM? 2+LWWb ) +1 . /■rbP22 + LWW'b"n ^ 

A60 = CK60*GO<^( 1 . / ( SM60 + LWOS ) +1 . 0 / ( SP60 + L WOS ) +1 . 0/ ( FRS + LWQS) ) 

— A ll S“s (_ < n 8 *GO* ( i . o / .1 bM I i 8 +L WOb ) + 1 . 0 / I bP 1 1 8 +-CWCT5T) ^ ^ 

A183 = C<183*GV»(l.o/(SM183+LWVS) + 1.0/(SP183-+LWVS) ) 

'\kt b = L<.Pb b*l)IVlZ*h W’^^LWW/ ( T Z ( 1 ) i . 6 ) 

AZ(I)= A183+A22+A60+A118+ARFS 

T71 coMTi'NijE ^ ^ irr: — : — 

DO R M=1.NM 

READ (5. 1 7) F.CM 7“ ^ ^ 

" it FORMAT ( 8F10.6 ) ' 

^ ^ ^ ^ 

KP=5 

prr=rn ^ — ; : 

DO PI 1=1 ,KR 

^ : — zLi=z( 1 r — : 1 ^ ^ 

IF (FR.FO.13.9) R( f) =PN*A1 39-»P( ZU)»*R139 

IF "T 'F R " ;FU. 8.9) FrfT7'HPN-»AByO.*PI ZU')*^^B8yO ^ ^ 

IF ( FR.FQ. 1 1 . 1 )■ R ( 1 ) =PN*A1 1 l»P( ZU)»*Blll 

p"! ~raTm xjttf — — — — ; — 

DO 1 1 I i<B .KT , . 

: iT^TZTrnTTi”; — :ir^: ~ ^ ^ 

IF( rU+8, ) .LT.0.01 ) GO TO 27 

■ CI^:=PN«FW**iVP5^fFXP r-6THb6« ( I' .+".'0TT43 «U )') ^ ^ 

GO TO 26 

2T-or^PrT^ R-«-TrT7-0‘(T7OrryPl“-H", ?/, l *'{ T; + IV1 C«~Li + ATr» 0 ^IIT) ^ ^ 

26 r M ) =rw>frK, . 
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n roNTiNi/p 
no 1 6 I -1 ♦<! 

wKi )♦ izrro'tKii ) .n i ) * a/a i ) : ,■ . ’■■ — 

? format ( 1 nx .'>( inx »Fi ?.a ) ) 

1 h A 7. f T ) = A z n I f T T^ ( f T > ^ ^ 

C THIS IS TO CALCULATE TR ANSM I T TANCE » TATM, TSK.Y 

UU 101 1H=1,1 ^ ^ 

H=H(IH) , 

JJ = 1U,*H ^ ^ ^ 

AT = 0.0 , . ’ • ' 

DU" 20 J = 1.JJ ^ r— ^ 

AT = AT+0.1*A7.( J) 

-ZD CT7NTTNTF ^ ^ ^ 

WRITE (6 *603) . 

■W3 hORVIA I ( 1HU.3X* /HDFV-.«b LS tUax , 1 '3H I RANSMl I I ANCE » iOX » i IH I / A IMOSRHER » 

T1 5X *5HT/S<Y/ ) 

i;U 30" < = 1*10 ^ ^ : 

THETA1=FL0AT(K-1 )*PI10 

IF(K.EO.IO) THETA1 = 8.95*PI10 ’ ^ ^ 

CES=C0S(THETA1 ) 

A T Z (K)=EXP 7 -i./r. !; sa" i ) ^ 

ATTN(K)=in.*ALOG10(ATZ(K) ) - ' 

sr~=" 0.0 — 

00 40 I =1 » JJ 

E7~S"TrrO ^ ^ ^ — - 

00 50 J=I»JJ 

S2 = S2+0.1#AZ(J) ^ ^ ^ ^ 

50 CONTINUE 

53=E'XP(-1 ./CES*S2 ) ^ ^ ^ " — 

si=si+o.i*s3*TZ( n*AZ( n 

■ZTD rUNI INDE ^ ^ 

TATM(K)=1./CES*S1 

!^4=0.0 ^ ^ ^ ^ ^ ^ 

DO 70 1=1.300 , : , 

S5=0'.O ' 7 

DO 80 I J=1 *1 

S5=SS+0.1*AZ( I J) " ^ ^ ^ ^ ^ ^ 

80 CONTINUE 

s6=EXP(-i ^ ^ ^ 

S4 = C4 + 0 . l*S6*TZn )«AZ( n 

/o CUNITMUL ^ ^ ^ 

TS<Y( < ) =1 ./CES*54 

X3 = C!C-1) *10 : 

WRITE(6*3) K3.ATZIK) »TATM(K) *TSKY(K) *ATTN(K) 

■3 F ORM AT'13'XT 15*'4{bX*E16.8)/) 

30 CONTINUE 

fHIS' IS' 10 CALCULA1L' EMI S. FREE .TM1E3 ^ 

WRITF(6*605 ) ■ 

^tT5 rTmrTATlTHITrr6XrTlRT)Fr,T7EE3^ ♦ 1 1 H E MM'l 'SSTVl IY.13X.11HT7REFLECTED 

1.15X.10HT/MFASUREO /> 

DO 44 X = ] » lO 
no 44 N=1 »2 
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I A I ^ I ^ » 

44 CONTINUE 


1 I c. 1 o » / 1 ) I _ 

71 FORviAK 10(‘5XF6.2)//) 


DO 93 1=1 »KR 


A4 l I J = A4 I 1 ) -It I 1 ) 

R( I )=0.0 

rD r 'n Trrr jK 

DO 15 I=KR»KT 


C( I )=0.0 
T~ CON T I NUE 
8 CONTINUE 





21 D{ I ♦J)=F(J+1 )-F( J) 


DO 22 I=2»6 
NI=7-T 

DO 22 J=1»NI 

?? D( r*J)=D( I-lO+l )-D( 1-l.J) 
C=7U 

P=Fr 1 )+ D( ] »1 )*7U 
DO 26 I=2»6 

C=C* ( ZLI-FLOAT ( I-l ) ) / FLOAT ( I ) 
26 P=P+C»D( I »1 ) 

RETURN 
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A”8. A compufer program for compufing the scattering characteristics 
of rain. 


■5 nJFNT 757TTWr~^H I H- T 5 FNG 

$ OPTION FORTRAN, FRCNT/15/ 

■5 FORTRAN NDFC< ' ^ ' 

C THE PROGRAM IS DESIGNED TO CALCULATE THE SCAT.CHARAT. OF PRECIPIT/ 

RPAL CAM" ■ 

DIMFNSION AMX { 100 ) ,PMX ( 100 ) , SX( 301 ) , SS ( 30 1 ) SH ( 301 ♦ 19 ) , SP ( 301 ) 

DIMENSION hCM 10) , rWI 10) »PR( 10) 

DIMENSION ANX(100),RNX(100) 

i.umh'lKX "ANX ,HNX ' ^ ^ 

COMPLEX R-,RI ,AMX »BMX 

PI -^^rrANn. i*4. ^ 

WRITE(6,111) 

— Ill " FORMA Tt T R l ,5 0 H T H E 5CAT TE1RT NG CHARAC TE RIS I I C5 ~ 0f PRECIPI l A l IONS TT 

READ(5,11 ) NI,ND»DX 

!1 ■fURMAr(2TS,l'i0.4) ^ ^ 

IF(DX.LF.n.O) GO TO A1 

NU=5^Xn ^ ' 

NA=NL*ND+1 

WRITKrb»T2I DX,NI,ND*NA ^ ^ ^ 

12 FORMAT { lH0,20HnASIC INTERVAL DX=, FIO , 4 , / /IX , 

1 2UHINTkX'RA I ED I rMbb"nNl =TT3 ♦ // IX , 20HCTIANGE Uh IN'IHVL ND= , 1 2 ,7 / 1 XT' 

2 20HT0TAL INT. PTS NA=»I5,//) 

REA0li>»2n Nl- ,NT ,NP,N!< 

21 F0RMAT(4I5) 

READ'I5T77) I F- 0 (J I » J= 1 ,NF ) 

READ(5,22) (TR( J) ♦ J=1 ,NT) 

KKAD1T,':T7 5 IPH IJ) , J=1 ,NP) " ^ 

22 FORMAT { 1 OFfl. 4 ) 

DO n I1=1»NF 

FR=F0(I1) ' 

lam=w,tfr 

DK=2.*PI /LAM 

— DPK=7ffH¥p-i /Dr*#3'«'2 ;o 

DO 8 12=1, NT 

tw^trtttt 

RI=R(TW,LAM) 

DO HTT=T,NP 

P=PR(I3t 

WRlIKrB'iTJ P,IW,hW ^ 

1 FORMAT ( 1H0,3H P= , F 1 0« 4 , lOX 4 3HT W= » FI 0. 4 , lOX , 3HFR= , F 10 . 4 ) 

DP = H 2 .'0 /TD K ttP'# ,7T ■ ) ^ ^ ^ 

L=1 

iTSTTrn ^ — 

Dxx=o.n . 

DU 16"'xr=l ,NI1 ^ ‘ ^ ^ 

DXX=DXX+DX ' 

DU 16~NN=1,NL " 

L = L + 1 

X = x4.fjxx ^ 

FXA=2.*nP<»FXP( -DP«X I*DXX 
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T 

C 


TALL ■ 
TF ( «iMX ) 
T g> — CON I' i NU r 
rr 


f UFFAH I X"» W I , AM7% F1MX ♦ SMX ♦ S.M5VNT 
4,1 ,74»74 


41 »24,?4 


7- 4 -- ^ X (L )=SM X * FX A 


Sr. (L) = SMS*FXA 

TF — I NK..I-C;< nr- GO lu 31 ^ 

DO ?6 K=1»NK 

T H E aP i /ib,* F1 : 0 A T ( K-l > ^ 

CALL COEFPT IN»A,MX »BMX ♦THE»PI ♦SPT»£TP) 


IF(SPT) 41»75,75 
75 rONTlMur 


75 

TTT 

51 

T7T 


IF(c,TP) 41»75*?5 
5H (L ♦< ) = ( <;PT + STP) »FXA 
TONT INUF 
CONT INIIF 


\ONT INUF 


51s5(ND»NI »DP.SX,NA) 

S?sS(Nr)»NI ♦DF,55»NA) ^ ^ 

5A=S2/5] ( 

wFin-ib,]5) si*b:^,s4 

15 FORMAT ( 1H0»20HFXTINCT. COEF BEX= ,E16.8//1X» 

r ?UH5CA 7T FR t NG COE F B5C=>E16»8//lX»20HALB ED O O F SCAIV 

IF (NK.c-Q.OI GO TO 32 

IVR I ! M b»bb I ^ : 

66 FORMATnOX»15HSCAT. ANGLES *.5X» lOHMAGNI TUDE ♦//» 

O i T 1 7 K,= 1^NK. 

bft=float(k-i )* in. 

DO 18 LN=2*r!A = 

5P(LN)=SH(LN*K) I 

1 R"CUN I INUh ^ 

S3=S(ND»NI ♦DP»SP»NA)/S2 

WRr i E(b.l4 T~ BE T *S3 

14 FORMATI?(10XF10.4 )//) 

17 CON T TNUF ^ 

3? FONT INUF 

-R"CnNI INDh ^ 


GO TO 43 


W= »hlU .A7 ) 


41 WR I 1 K I 6,471 

47 FORMAT I 1H0*38HINC0RRECT INPUT OR CALCULATED RESULTS ») 

43 Slop 

END 

FTTTrn770T~NT7FrF 

THE FUNCTION INTEGRATES SCAT . CHARACT . BY COMPOSITE RULES 

FUNC I I ON S f Nl.' » N r , OP'. PZl N A ) ^ 

dimension F( 6 )»PZ(NA) 



E=n,n 

r ^3A= I . .n ^ 

F (6 ) =n,o 

DO- 17 I. L = I rATD : 

DO 73 n =1 ,NI 

FnT^TT 6 -l 

DO 1 6 L =2 *6 ( 
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16 F(L)=PZ{J) 

ARC=FL0AT(LL+1 ) /FLOAT (LL) 

/ arc 

17 F(6)=ARr»F(6) 
c,= c,'+f (6T*("RA/r>P 
RPTURM 

“FTTP 

% FORTRAN NDFr< 

t ,IHE SUBROUl irjE CAL'CULAFtb 1 Ht M 1 E» EX T I NC T . AND' SCA T . CO'EV 

SUBROIITTNF C0FFAP(X»RI ,AMX,RMX»SMX,6MS*N) 

COMPLEX Yo R1 *A'/';'AVU,WXNi,WXU,W'X*AAX»BAX,AMX'7BP0<,CAX,bDX 

DIMENSION AY{ 100 ) »WX ( 1 00 ) *AMX ( 100 ) »BMX ( 100 I 

yr>r«-Ri 

YR=RFAL(Y) 

YI=-AIMAG(Y) ^ ” 

YNI=.SINfYR)«-rOSfYR) 

YL=Z.»Y I 

YN2=(FXP(YL)-FXP(-YLn*.2 5 
' YL) = SIN(YK)*51N(YK)-o,b + 0.25*( EXP( YL ) +tXP ( -YL ) > 

AY0 = { YN1+ (0. .1. )*YN2 ) /YD 
AY ( 1 ) = -l ,/Y + l . / I 1 ,"/Y-AY0 ) 

V'XO = SI N< X ) + ( 0. »1 . ) »COS{ X ) 

WXN]=CUS(X)-(0.»1. )*SlN(X) ' 

WXIl )=WX0/X-WXN1 
WXi 2 ) =T,»WX ( 1 )/X-WXn 
CAX = AY(1 )/RI+l./X 

pox'll *AYn.")+i .XX 

AMX ( 1 ) = (CAX*RFAL( WX( 1 ) )-RFAL ( WXO ) ) /(CAX*WX( 1 )-WX0) 

PMXI 1 ) = (BDX*RFALI'C^fl ) ) -RFAU W:»<0') ) '/ ( BDT*WX ( 1 )'=W7D } 

ACR=1 .n'^-06*CAR,S ( AMX ( 1 ) ) 

RCI^=l.nF-06*CAPS(RMX( 1 ) ) 

IF(ACR.LT.1.0E-06) ACR=1.0F-06 
IF ( BCR.LT. 1 .OE-06 )'”B(:R=r. 0 E -06 
S'<X = 0.0 

SMS = 0.0 

DO 3 1=1,50 

IF ( I . fo. 1 ) GO TO ^ ^ ; : 

IF (I.F0.2) GO TO ? 

A I =2* I -1 ; 

WX( I ) = AT*WXf I-l')/X-WX(I-?) 

7 — 3T5-I 

AY(I)= -AI/Y + 1./(AI/Y-AY( I-l ) ) 

AAX = AY( 1 I / RI 4-A I /X 

RAX = RI*AY( I l+M/X 

AMXn ) = I AAX>RtAL ( WX ( 1 ) I -REAL ( WX( 1-1 ) I ) / ( AA■X^^ WXTl T-WX (1-11) 

RMXd ) = (RAX»REAL(WX( 1 ) )-REAL(WX( I-l ) ) )/(3AX*WX( I ) -WX ( I-l) ) 
IF(RFAL(AMX( 1 ) +bM a ( I ) ).LE.1.L-1U) GO I 0~6 ^ 

9 continue: 

zrr^irm 

c;my = snX + AI*RFAL( AMX ( I ) +BMX ( I ) ) 

«TMS = SMS + AT* ( GAMS ( A’MmT) «-*2 + rAR'SmM'Xl I ) ) 

AMA = CARS ( AMX ( 1) ) 

MMA'^ATISIRMX (III 

IF(AMA.IF.ACR.AND.BMA.lE.RCR) go to 6 
3~rnTrrrwnF 


142 




— - — 

% FORTRAN NOFCK 

r T H F SUB R OU TI N E CALCUL A TE'S 1 HF THAR AC.'CURVFS "UF bC. A ( « Vb. ANbCtS" ' 

SURROUTINF C0FFPT(N»AMX»RMX»THE»P1 »SPT *STP) 

CDMP I.FX AMX ^RMXtr PT ^f, TP ^ 

OIMFNSION AMX( 1 00 ) »^MX (1 00 ) ♦fPT ( 100 ) »CTP ( 1 00 ) »PH I (1 00 ) »TAU ( 1 00 ) 

phi: 1 ) = r.o “ ^ ; 

PHI t?)=0.0»rOS(THF) 

P A t !)n I ^T DS IT H F ) 

TAi;:2 »=?.0*C05(2.*THF) 

e ;N 2 = F T NC T HFr* ;= ; T Nf T H F ) — 

r=T=(o.o»o.o) 

C T P=( 0 .0 >0.0 1 ^ 

DO 3 1=1 >N 

:: TF T i .r: 'p~ .2 I i>o U) 9 ' — " — 

AT=2*I-1 

Bi = r - i 

PHI ( T ) = (TAUn )*AT*PHl ( I-l l-PHI ( 1-2) )/BI 

TAUn ) =TAU( 1 ) *T PH n I ) -PHI (1-2) )-AI*SN2*PH I ( I - 1 ) *T AU (I -2 > 

9 CONTINUE 

Cl =FL0AT (2*1 + 1 ) /FLOAT ( I*( I+l ) ) 

CPT = CPT+CI*( AMX( I )*PHI ( I ) +R(«1X ( I )*TAU( I ) ) 

Ctk = CTP + CI*CaMX( I )«tAU( I )+BMX( I )*PHl'n ) ) 

3 CONTINUF 

SPT=CARS (CPT ) *CARS (CPT ) / (A.*PI ) 

STP = CARS (CTP)*rARS(CTP ) /(A.*PI ) 

RETURN 

FND 

3K PflRlRAM NDFrt< 

C THE FUNCTION CALCULATES THE INDEX OF REFRACTION OF WATER 

FUNCl I ON rttr ;r'am ) — ^ 

COMPLEX F>R 

Re ftL R AM »h0 >E INF >DLAM "" 

FINE = 3.5 

DT=7'R-273. 

EO = 87. 7A-0.A008«DT+9. 398*1 .E-A*DT*DT+1. 41*1. E-6*DT*DT*DT 

OLAM=l. ilB+O'; 1 122#L) I + i .8*r.'E-^7^Ur*UT-8". 3*r.L-6'#Dr*DT*0T 

F= (EO-F INF )* ( 1 . >0. ) / ( ( 1 . >0. )+(0. >1 . ) * (DLAM/KAM) )+EINF 

R=rSDRT (I- ) 

RPTURN 

Fxrn 

T. FXECUTF 

LTMi Ts. — rrr > 20000 > >"i ^>000 
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A-9. 


A compufer program for computing the apparent temperature 
of heavy rain . 


MDENl 7576TWTT7TTFrm ^ 

f. OPTION FORTRAN, rpcNT/I 3/ 

% FORTRAN Mr>rf-|^ ^ 

OIMFNSION A(3,5),B(3,5),C(5,5)»0(5<.5)»P(5»1),C0(3»3)»P0(5»31» 
ICA ( 3»3 ) ,PF { 3 , 3 ) ,RA ( 3»3 I TrrnTF , 3T';X Hj') ( 5'f , L ( t ,'M( 3 > ,2 ( 3 » 3'uT, 

?RO (3,3)»RF(3,3),RF{3»5)»Pr(3»l)»RC(3»3)»CC(3»5)»CR(3,3)»C3{3»3)* 

^ •5PPt3Tl ) ,r;Tt 3 , M »FP ( 3M ) »T'D( 5 , I riOrr 3 ,'5 , AO ) ,PHLU( 3 ) » bL( 3»3,"iO IT 

4DL (3, 30), PL (3, 30) ,PHP ( 3 , 30 ) , PHPO ( 3 ) » PLO ( 3 ) , QLO ( 3 , 3 ) , PHL ( 5 » 3 C ) , 

5 ^ PT5 ,T) ,HPl 5 , 1 r »FB( 3, 1 ) »(_P ( b'Vl )TFC r'3, n »PA{ 3 ♦ 1) »P!3 ( :> , IT 

COMMON U(5 ) ,W(5) 

READ (5, 31) U,W ^ 

51 FORMAT ( 5F10. 3 ) 

■ WRlTb(6,T) UVW ^ 

1 FORMAT ( 3 ( 3XF16.B )// ) 

~Dir IT LL"=i i7 ^ : ^ 

REAO(3,5) ARD»P0,P1 ,P2,TAU,TW,TG 
b' FO'RMAT (7TlriT3 ) ^ ^ ^ 

WRITE(6,3) ABn»P0,Pl »P2»TAU,TW,TG 

CALL AHCI- ( AHD,P0,P1 ,P2 ♦ I AU, I W,A»H »(. »h'i 1 ) 

WRITE(6,1) A,R,C»F 
^ UO 1 3 1 L = 1 ♦ 4 ^ 

READ(3,5?) Q«P 

5'2 fORMAT(3T10'.5 ) 

WRITE(6,1) Q»P 

DO 11 1 = 1,5 

P{ I ,1 ) =P ( I , 1 ) «-TG 

— puum — =Pf~r,'iT ^ ^ ^ 

DO 11 J=l,3 

OLO( 1 , J ) =0( 1 , J ) 

11 CONTINUF 

^ ^ ^ 

2 CALL MPRD(R, 0,30,3,3, 5) 

CTTO:: M ATTD ( TIO ', CV H Cl 3 

CALL MPRD(C,0,CQ,3,3 ,5 ) 

CATE MADD(CO,A,CA, 5,»51 : 

CALL MINV(CA,3,D,L,Ml 

CATX MPRD(RC,CA,R'A,3,5,5) ^ 

CALL MPRDIRA, R,RB,3,3,5) 

( ALL 7'1'bimTRfTV CTRC; 3 jE'I ^ 

CALL MPRn(RA,C,Rn,3,3,3) 

CALL 'MS MET 7TJ RDTRF",'3 ♦ET 

CALL MINV(RE,5,D,L,M> 

CALL MPRD ( RF; R(^ R FI 3 »'5 751 

CALL MPRD(B,P,BP,3,3»1 ) 

C AL L w s I ! |3 ( L , R P 7F 3 ,'5",' H ^ 

CALL MPRD(C,P,rP,5,5»l ) 

CALL M ,3 { IB'I FTC P 7FC ,“5“,' 11 

CALL MPRD ( RA ,FC ,PA , 3* 3 , 1 ) 

^ CALL F G 1 1 ■( 1 H' ; P A", P RT5 » T3 ^ ^ ' 

CALL MPRDIRF ,PR,PC ,3»3 , 1 ) 

^ CALL MPR[)( rVRF,( R ,3* 3',5T~^ ' 

CALL MADOtCR, R,rR,3»3) 
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C ALL "> ^ PR n r rA »cn >rr , 5*3 <, 5 ) 
CALL VAOnCPr, P»PP.5»1) 
C ALL M PR n t — r ,PP ,CT *5»5 HT 
CALL MSUB( F»CT*FP,5*1) 
CALL MPR!’RA»TT’*L)I>»5»'5, 1 ) 
DO A I =1 »5 

DL ( T » r J) ^DDf r .i) 

PL ( I * IJ) =PC ( I . 1 ) 


DO 9 J=1 »5 

OL ( I »J ♦ IJ ) =RF ( I ) 

GL ( I ♦ J» IJ )=CC ( I » J ) «(-l , ) 

0( I , J) =RF ( I ,J ) 

9 CONTINUE 

P( I »1 ) =PC( 1*1) 

“4 CONTIMuf 

IJ=IJ+1 

TF n J-'ioT~2T2" ^ 

3 CONTINUE 

Du b 1=1*3 ^ 

PHLII* 30)=0.0 

DO 7 JJ = 1*3U 

K=31-JJ 

IFK.KU. IT bU 10 12 

21 T,iO=n.O 

DO' 31 J= 1 T 3 

31 Z^( U<)=2( I,K)+GL( I »J*K:)*PHL(J*ia 
Pm LI 1 ♦><; - 1 ) rZ ( F ,<) +f)L ( 1 * < ) 

12 CONTINUE 

ZmXT=TTT75 ^ 

DO 32 J=l*5 

/:u*<)=7(I,<) +0L ( I ,J*< ) *PHL ( j*k. ) 

PHP( I ♦< ) =Z ( I *<)+PL (1*0 

■ 7 CONI INUb 

X( I )=0.0 
Du 3 3 J = 1 * 3 

33 X( n=X( 1 )+GL( 1 *J*1 )*PHL( J*1 ) 

— pHLo n ) =xTiT+T)TTrrr? ^ 

Y( I )=o.n 

L'U'”3A ■ J=1 73 

^4 Y( n=Y( I )+0L0( I *J)»PHL0( J) 

PhPU( I ) =y I ! )+HLU( I ) ^ 

6 CONTINUE 

W RITE (6,1 )' P H P 0 ' * P H P * P HL0 — " ,'PHL 

13 CONTINUE 

STOP 

FND 

FORTR AN ND F CK 

GURROUTINF ARCF ( ARD * PO * P I * P2 * T AU * TW ♦ A , P ,C * F , K ) 

DT ^ F M G T OAT A ( 5 ♦ 3 )- VPT ' 3 , 5 ) « C T 3T 3 ) F 

COMMON U(3),W(3) 

-D O ' " ? T = 1 * 3 

IF(K. r.F.?) TW = 0.0 

F f T * 1 ) = n r-ARm > t w 

DO 2 J=1 *3 

ro T = - . 3-»^~A-prr)Krf7n 

GO TO n 1 .22,33) ,K 
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T 


' /• i i ( I. M 1 / — I . » ; 


PJ-'. ( ■^ ,*U (J I »IJ( J ) -1 


rOF*(Pn + Pl>HM I ) »U( J)+P?*PI^^PJ) 


GO TO 


I \ L • - U l I ) 


PJ1=50RT (1 .-(J{ J)«ur J) ) 


PJ2=3.*U(J)*PJ1 


I L 9 sJ f 


B( I »J)=C0F*(P1*P11«PJ1+P2/3.*PI2*PJ2) 



RETURN! 

END 


EljP^PljTlNE MSU?( 3, 

D I ^ E • ' S I n V A f N , y ) , R ( N , M ) I R ( M , M ) 
I" D 1 I = i » N 
PO 1 J=1»M 

1 P ( ! . J ) = A ( I , J ) - T { I , J ) 

Ft; turn 
END 


SU-K'^UTINE ;UUP( A, 3,R|f>)*H) 
r I N'f: MS I aM A ( 'I . R ( N , M ) , R ( \l , y ) 
PQ 1 1=1. N ■ 

ro 1 j=i.M 

^'pTljRlj 

LImT 

EU-’lvOUTlME iMP!n(A,3,R,f.!,M,L) 
n I "l/'S I ON A ( M , y ),[•)( M , L ), R ( M , L ) 
P'O 1 I = 1.N . 

P() 1 t; = l,L 
F (! . K ) = ;) , 0 
r* 0 1 J = I , M 

I "R (-r , K) = ;^( j ,K ) 4 -\ ( 1 , J )»Q ( j, K ) 
RRTliRI. 

rsr 

* EXEGUTF 

T. LIMITS 10.30000. .1.^000 
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